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’  yg  report  experimental  response  functions  at  twelve  gamma- 

ray  energies  over  the  range  of  123.6  keV  to  11.67  MeV  for  a  2 . 54-cm-diameter  x  2.54- 
cm-long  and  a  7 .62-cm-diameter  x  7.62-cm  long  bismuth  germanate  scintillation  detec¬ 
tor.  The  measurements  were  made  at,  or  corrected  to  correspond  to,  source-to-detec- 
tor  distances  of  100  cm.  The  resolutions  of  the  detectors  at  a  gamma-ray  energy  of 
661.6  keV  are  11.6  percent  for  the  small  detector  and  13.0  percent  for  the  large 
detector .  ^4Je  also  present  an  interpolation  method  for  generating  a  response  function 
at  any  gamJJfe-ray  energy  in  the  range  of  123.6  keV  to  11.67  MeV  for  either  scintil¬ 
lator  using  the  experimental  response  functions.  Additionally,  this  method  is  used 
in  constructing  response  matrices  for  unfolding  gamma-ray  pulse-height  distributions 
acquired  with  the  detectors.  A  computer  code  written  in  FORTRAN-77  generates  the 
response  matrices. 
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INTRODUCTION 


I . 

As  part  of  an  extensive  experimental  study  of  the 
energy  and  angular  distributions  of  neutrons  and  gamma  rays 
produced  from  50-  to  200-MeV  proton  and  deuteron  bombardment 
of  thick  elemental  targets,  we  acguired  two  sizes  of  bismuth 
germanate  (BGO)  inorganic  scintillators,  2 . 54-cm-diameter  >; 
2.54-cm-long  and  7 . 62-cm-diameter  x  7 . 62-cm-long ,  for  chs 
measurement  of  gamma-ray  energy  spectra.  The  gamma-ray 
energies  of  interest  ranged  from  approximately  250  keV  to  15 
MeV  with  the  most  intense  part  of  each  gamma-ray  spectrum 
lying  below  E-y  =  4500  keV  (Refs.  1-3).  The  neutron  energy 
spectra  are  dominated  by  neutrons  with  Ej^  <  10  MeV  (Refs.  2- 
6)  . 

BGO  was  chosen  instead  of  Nal(Tl)  because  of  its 
favorable  characteristics  regarding  detection  of  high-energy 
gamma  rays  in  the  presence  of  both  large  numbers  of  low- 
energy  gamma  rays  and  large  numbers  of  neutrons  with 
energies  primarily  less  than  10  MeV.  In  one  of  the  early 
investigations  of  the  gamma-ray  response  of  a  BGO  detector, 
Evans  (Ref.  7)  found  the  ratio  of  the  photopeak  efficiency 
of  a  3.8-cm  X  3.8-cm  BGO  scintillator  to  that  of  the  same 
size  Nal(Tl)  scintillator  to  be  3.3,  4.5,  and  5.6  at  E-y  = 
662,  1332,  and  2754  keV,  respectively.  He  concludes  that 
because  of  the  increased  photopeak  efficiency  and  the 
relative  flatness  of  the  BGO  photopeak  efficiency  curve  as 
compared  to  Nal(Tl),  BGO  detectors  are  the  best  choice  for 
measurement  of  high-energy  gamma  rays  accompanied  by  large 
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numbers  of  low-energy  gamma  rays.  Further  studies  of  BGO 
detectors  (Refs.  8-10)  show  that  the  neutron  sensitivity  of 
a  7.62-cm  X  7.62-cm  BGO  scintillator  is  considerably  less 
than  that  of  the  same  size  Nal(Tl)  detector.  Hausser,  et  al 
(Ref.  9)  show  that  the  BGO  to  Nal(Tl)  neutron  sensitivity 
ratio  is  approximately  0.2  at  Ej^  =  500  keV,  1  at  Ej^  =  3500 
keV,  and  1.2  at  E^^  >  5000  keV.  The  authors  of  Ref.  9 
predict  that  this  ratio  is  at  least  less  than  0  1  for  Ej^  < 
500  keV  and  conclude  that  BGO  as  compared  to  Nal(Tl)  has  a 
much  greater  gamma-ray-to-neutron  detection  ratio.  The  only 
disadvantage  of  a  BGO  scintillator,  as  compared  to  Nal(Tl) 
discussed  in  the  studies,  is  the  poorer  energy  resolution  of 
BGO  (Refs.  1  and  4).  However,  in  our  work,  energy 
resolution  of  the  gamma-ray  detector  is  not  important  since 
the  gamma  rays  produced  from  medium-energy  charged  particle 
bombardment  of  targets  primarily  form  a  continuum  (Refs.  2 
and  3 ) . 

In  order  to  obtain  gamma-ray  energy  spectra  using  a  BGO 
detector,  the  detector  response  must  be  deconvoluted  from 
the  measured  gamma-ray  pulse-height  distributions.  The 
deconvolution  process  is  commonly  performed  with  response 
matrix  unfolding  computer  codes  such  as  the  widely  used 
FERDO  code  (Ref.  11).  The  response  matrix  input  to  the 
unfolding  code  consists  of  numerous  BGO  response  functions 
spanning  the  gamma-ray  energy  range  of  interest.  A  response 
function  is  the  pulse-height  distribution  resulting  from  a 
monoenergetic  gamma-ray  source  of  unit  intensity  incident  on 
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the  BGO  scintillator.  Each  response  function  can  be  quite 
complicated  in  form  and  depends  on  many  factors;  the 
quality  of  the  BGO  scintillation  crystal,  the  photomulti¬ 
plier  tube  and  associated  detector  electronics,  the  gamma- 
ray  source  and  detector  configuration,  and  the  detector 
shielding. 

The  standard  practice  for  obtaining  a  set  of  response 
functions  is  to  calculate  or  measure  them  for  a  BGO  detector 
in  a  specific  experimental  arrangement.  The  calculations 
are  performed  with  Monte  Carlo  electron-photon  transport 
codes  such  as  CYLTRAN  (Ref.  12)  and  PEP-C  (Ref.  13). 

Although  these  codes  produce  photopeak  efficiencies  which 
agree  well  with  measured  efficiencies  for  different  sized 
BGO  detectors  (Refs.  13  and  14),  there  is  no  extensive 
comparison  between  calculated  and  measured  response 
functions.  Numerous  measurements  of  photopeak  efficiencies, 
total  efficiencies,  resolutions,  multigamma  responses,  and 
response  functions  for  2.54-cm  x  2.54-cm  and  7.62-cm  x 
7.62-cm  BGO  detectors  exist  (Refs.  8  and  10,  13-17). 

However,  the  only  complete  set  of  response  functions  for  a 
single  BGO  detector  and  gamma-ray  source  configuration  is 
the  result  of  the  measurements  of  Moss,  et  al  (Refs.  16  and 
17)  for  a  7.62-cm  x  7.62-cm  detector  positioned  30  cm  from 
various  gamma-ray  point  sources.  These  response  functions 
are  17-parameter  analytical  functions  fitted  to  some  of  the 
measured  pulse-height  distributions.  For  our  work  we  found 
tnese  theoretical  response  functions  to  be  inadequate. 
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especially  at  the  higher  gamma-ray  energies,  because  they  do 
not  include  several  effects:  (1)  the  energy  dependence  of 
the  detector  resolution  within  a  given  response,  (2)  the 
asymmetries  of  the  photopeak,  single-escape  peak,  and 
double-escape  peak  due  to  secondary  radiation  and  charged 
particle  escape  from  the  crystal,  (3)  the  variation  in  shape 
of  the  Compton  distribution  with  increasing  gamma-ray 
energy,  and  (4)  in-scattering  of  gamma-rays  from  the 
detector  housing  materials  and  surrounding  shielding 
materials.  These  effects  vary  in  importance  depending  on 
the  BGO  detector  and  the  experimental  arrangement.  We  were 
primarily  concerned  with  the  first  three. 

Because  of  the  inadequacies  of  the  calculated  and 
measured  response  functions  of  BGO  scintillators  reported  in 
the  literature,  we  measured  the  response  functions  for  a 
2.54-cra  X  2.54-cm  and  a  7.62-cm  x  7.62-cra  BGO  detector  at 
twelve  gamma-ray  energies  between  123.6  keV  and  11.67  MeV. 
This  paper  reports  these  measurements  along  with  tabulations 
and  graphs  of  intrinsic  photopeak  efficiencies  and  detector 
energy  resolutions.  Furthermore,  we  present  a  computer  code 
written  in  FORTRAN-77  which  uses  these  measured  response 
functions  to  generate  a  response  matrix  for  use  in  unfolding 
codes.  The  resultant  response  matrix  can  be  used  for 
unfolding  gamma-ray  pulse-height  distributions  over  the 
energy  range  of  approximately  600  keV  to  12  MeV.  Both  the 
reported  response  functions  and  response  matrix  generation 
code  should  be  of  general  use  because  they  can  be  used  for 
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other  BGO  detectors  of  the  same  size  with  equivalent  or 
poorer  energy  resolution  merely  by  "smearing"  these  response 
functions  with  the  appropriate  resolution  factor. 
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II.  RESPONSE  FUNCTION  MEASUREMENT  DETAILS 


BGO  Detectors 


We  used  2.54-cin  x  2.54-citi  and  7.62-cin  x  7.62-cin  BGO 
detectors  manufactured  by  Bicron  Corporation.  Each  detector 
came  as  an  integral  assembly  with  photomultiplier  tube  and 
magnetic/light  shield.  A  schematic  diagram  of  each  detector 
is  shown  in  Fig.  1  and  Fig.  2.  Not  all  of  the  packaging 
details  of  the  scintillators  are  shown  in  the  figures. 
Further  details  can  be  obtained  from  Bicron  Corporation  by 
specifying  the  model  numbers:  2.54-cm  x  2.54-cm,  Model  # 
lMlBGO/1.5;  7.62-cm  x  7.62-cm,  Model  #  3M3BGO/3. 
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Figure  1.  schematic  of  the  2 . 54-cm-diameter  x  2.54-cm-long 
BGO  scintillation  detector.  (Schematic  shown  with  permis- 
ion  of  Bicron  Corporation). 
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Figure  2.  Schematic  of  the  7 . 62-cm-diameter  x  7.62-cm-long 
BGO  scintillation  detector.  (Schematic  shown  with  permis¬ 
sion  of  Bicron  Corporation). 


Gamma-Ray  Sources 

For  the  response  function  measurements  we  used  eight 
standard  gamma-ray  point  sources  of  known  intensity  spanning 
the  energy  range  of  123.6  to  1836  keV  (Table  1).  We  also 
used  four  sources  of  gamma  rays  of  unknown  intensity  from 
the  ^^Na(n,'y)2^Na,  ^Be  («  ,  n -y )  ,  ^^F(  p,  or-y  )  ,  and 

^^B(p,'-y)^2r  reactions  spanning  the  energy  range  of  2754  keV 
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Table  1.  Gamma-Ray  Sources  Used  For  BGO  Scintillation 
Detector  Response  Function  Measurements. 


Nuclide  or 

Reaction 

Gamma-Ray  Energies 
(keV) 

57co  (a) 

122,  136.5  }123.6 

203Hg  (a) 

279.2 

85si-  (a;' 

514 

137cs  (a) 

661.6 

54Mn  (a) 

834.8 

65zn  (a) 

511,  1115.5 

^^Na 

511,  1274.5 

88  Y  ( S' ) 

898,  1836 

23Na(n,'y  )24Na  (t)) 

1368.6,  2754 

^Beicx,n'r)^^C 

others  <  1000,  2223,  4440 

1^F(P,c»:'y)^^0 

511,  6140,  6920,  7120 

11b(p,t)12c  «J) 

511,  4440,  11,670 

Calibrated  point  source  purchased  from  Isotope  Products 
Laboratory . 

Uncalibrated  extended  source  produced  by  neutron 
activation.  ♦ 

Uncalibrated  extended  ^^^Pu/Be  neutron/gamma-ray  source. 
Uncalibrated  expended  source  produced  by  radiative 
proton  capture . 


to  11.67  MeV  (Table  1).  The  point  sources  had  a  nominal 
strength  of  100  >L/.Ci  and  were  of  a  thin,  "scatterless”  disc- 
type  construction  (Ref.  18).  The  uncertainty  in  intensity 
for  each  source  was  certified  to  at  most  ±  4.6  percent. 

Each  of  the  gamma-ray  sources  produced  from  neutron,  alpha, 
or  proton  reactions  was  an  extended  source  (Table  2).  The 
reaction  ^^Na(n,'y)2^Na  was  produced  using  sodium  motal 
coated  with  paraffin  and  bombarded  with  thermal  neutrons 
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Table  2.  Summary  Of  Extended  Gamma-Ray  Sources  Dimensions. 


Nuclide  Designation 
and  Reaction 

Source  Description  and 
Dimensions 

^^Na(n,'y  )^^Na 

2.5-cm  X  2.5-cm  square,  1  cm 
thick 

238pu/Be 

^Be(Q:,n'Y)^2c 

0 . 8-cm-diameter  x  3-cm-long 
cylinder  in  large  water  tank 

16o* 

^^F(p,ca:'y)  ^^0 

2-cm-diameter  x  0.4-cm-thick 
disk 

12c* 

11b(p,-,^)12c 

1 . 2-cm-diameter  x  0.2-cm-thick 
disk 

The  large  water  tank  served  as  a  neutron  moderator. 


from  the  United  States  Air  Force  Academy  (USAFA)  200-keV 
Cockroft-Walton  accelerator /neutron  generator.  Both  the 
^^F(p,c»:'y ) and  reactions  were  produced  using 

the  USAFA  400-keV  Van  de  Graaff  accelerator  with  proton 
energies  somewhat  greater  than  the  resonance  energies  at  Ep 
=  340  keV  and  163  keV,  respectively.  The  targets  consisted 
of  pressure-packed  powder  discs  of  sodium  fluoride  and 
boron . 

Experimental  Arrangement 

For  the  response  function  measurements,  the  eight  point 
sources  were  on  the  detector  axis  line  a  distance  of  100  cm 
from  the  front  face  of  each  BGO  scintillator.  The  choice  of 
this  distance  was  a  compromise  between  the  data  acquisition 
rate  driven  by  source  intensity  and  the  desire  to  have  a 
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nearly  parallel  beam  of  gamma  rays  incident  on  the  detector. 
The  center  of  the  cylindrical  ^^^Pu/Be  source  was  also  100 
cm  from  the  front  face  of  the  BGO  detector.  However, 
between  the  source  and  detector,  10  cm  from  the  end  of  the 
238pu/Be  cylinder,  we  placed  a  3 . 81-cm-diameter  cylindrical 
plug  consisting  of  a  3-cm- thick  piece  of  polyethylene,  a 
0.9-cm- thick  piece  of  lead,  and  a  4-cm-thick  piece  of 
polyethylene  (Fig.  3).  This  plug  reduced  the  number  of 
neutrons  and  low-energy  gamma-rays  from  the  source.  The 


Figure  3.  Experimental  arrangement  for  the  BGO  response 
function  measurements  using  the  2^®Pu/Be  neutron/gamma-ray 
source . 
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outside  of  the  water  tank  facing  the  detector  was  lined  with 
5.08  cm  of  lead  to  reduce  the  number  of  2223-keV  capture 
gamma  rays  escaping  from  the  tank.  The  center  of  the  ^^Na 
source  was  located  5  cm  from  the  front  face  of  the  2.54-cm  x 
2.54-cm  BGO  detector  and  10  cm  from  the  front  face  of  the 
7.62-cm  X  7.62-cm  detector  due  to  the  extremely  low 
intensity  of  the  ^'^Na  source.  The  centers  of  the  sodium 
fluoride  and  boron  source  targets  were  34  cm  from  the  front 
faces  of  the  detectors.  The  targets  were  positioned  at  the 
middle  of  an  evacuated  52 . 1-cm-diameter  aluminum  scattering 
chamber  with  thin.  Plexiglass  viewing  ports  (Fig.  4). 


Figure  4.  Experimental  arrangement  for  the  BGO  response 
function  measurements  using  the  sodium  fluoride  and  boron 
source  targets. 
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Again,  this  distance  was  chosen  due  to  the  low  intensity  of 
the  sources.  Since  the  response  function  for  a 
scintillation  detector  is  known  to  vary  with  source-to- 
detector  distance  (Ref.  19),  the  point  source  was  also 

used  for  response  function  measurements  at  distances  of  5  cm 
and  25  cm  for  the  2.54-cm  x  2.54-cm  BGO  detector  and  10  cm 
and  50  cm  for  the  7.62-cm  x  7.62-cm  BGO  detector.  These 
measurements  along  with  further  information  from  the 
literature  was  used  to  modify  those  measurements  taken  with 
source-to-detector  distances  less  than  100  cm  in  order  to 
make  them  correspond  to  100-cm  source-to-detector  distances. 
This  modification  is  described  in  the  analysis  section  of 
the  paper.  A  summary  of  the  source-to-detector  distances 
for  the  response  function  measurements  is  presented  in  Table 
3. 

Gamma-ray  scattering  effects  were  minimized  for  as  many 
of  the  response  function  measurements  as  possible.  Both  the 
detectors  and  sources  were  a  minimum  of  100  cm  from  the 
floor,  ceiling,  and  nearest  wall.  The  eight  manufactured 
sources  along  with  the  ^^Na  source  were  suspended  with  tape 
from  the  ceiling,  hanging  150  cm  above  the  floor.  Both  BGO 
detectors  were  mounted  on  aluminum  camera  tripods  150  cm 
above  the  floor.  The  scattering  effects  were  larger  for  the 
remaining  sources  since  the  238p^/5g  source  was  contained  in 
a  water  tank  and  the  accelerator-produced  sources  were 
targets  mounted  in  an  aluminum  scattering  chamber.  However, 
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Table  3.  Summary  Of  Gamma-Ray  Source-To-Detector  Distances 
For  The  2.54-cm  x  2.54-cm  And  7.62-cm  x  7.62-cm  BGO 
Scintillation  Detectors. 


Nuclide  or 
Reaction 

Source-to- 

Detector 

Distance 

2.54-cm  X  2.54-cm 

(cm) 

Source-to- 

Detector 

Distance 

7.62-cm  X  7.62-cm 

(cm) 

0 

o 

IT) 

100 

100 

203Hg 

100 

100 

85sr 

100 

100 

100 

100 

54Mn 

100 

100 

®^Zn 

100 

100 

22Na 

100 

100 

00 

00 

>< 

5,  25,  100 

10,  50,  100 

^^Na(n,-y’)^^Na 

5 

10 

^Be  ( cK ,  n-y )  ^ 

100 

100 

^9F(p,Q5'r)i^O 

34 

33.7 

llB(p,'y)^2c 

34 

33.7 

All  distances  are  measured  from  the  center  of  the  gamma- 
ray  source  to  the  front  face  of  the  BGO  scintillation 
detector. 


the  detectors  remained  mounted  on  aluminum  camera  tripods 
for  these  measurements. 

We  simultaneously  acquired  pulse-height  distributions 
for  both  the  2.54-cm  x  2.54-cm  and  7.62-cm  x  7.62-cm  BGO 
detectors  for  each  gamma-ray  source.  The  dynode  signals 
from  the  photomultiplier  tube  bases  were  fed  through  ORTEC 
113  preamplifiers  into  ORTEC  450  or  572  amplifiers.  The 
resulting  amplified  signals  were  digitized  into  1024 
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channels  with  Nuclear  Data  65  and  66  multichannel  analyzers. 
For  each  pulse-height  distribution  measurement,  we  adjusted 
the  amplifier  gains  to  maximize  the  1024-channel 
digitization.  The  pulse-height  distribution  acquisition 
times  ranged  from  1  to  5.5  hours.  Before  and  after  each 
measurement,  we  performed  an  energy  calibration  of  the 
pulse-height  axis  using  4  to  8  known  gamma-ray  photopeaks. 
There  were  no  noticeable  gain  drifts  during  any  of  the  data 
runs.  As  soon  as  possible  after  a  data  run,  we  acquired 
background  pulse-height  distributions  by  removing  the  source 
from  the  laboratory  room.  For  most  of  the  measurements,  the 
background  run  times  were  equal  to  or  longer  than  the  data 
run  times.  Furthermore,  the  pulse-height  axes  of  the 
background  distributions  were  energy  calibrated  separately 
from  the  data  runs.  Both  the  data  and  background  1024- 
channel  distributions  along  with  the  energy  calibration 
spectra  were  stored  in  ASCII  format  for  analysis  using  a 
personal  computer  and  a  MicroVAX. 
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III.  ANALYSIS 


The  first  step  in  converting  the  measured  pulse-height 
distributions  into  response  functions  was  energy  calibration 
of  the  pulse-height  axis  of  each  data  and  background  file. 

We  fitted  each  peak  of  a  calibration  spectrum  with  a 
Gaussian  function  plus  a  second-order  polynomial  function 
using  the  generalized  function  minimization  program  MINUIT 
(Ref.  20).  For  closely-spaced  peaks,  simultaneous  fits 
involving  multiple  Gaussian  functions  plus  one  second-order 
polynomial  function  were  used.  The  result  of  these  fits 
were  specific  energy  values  corresponding  to  channel 
numbers.  We  then  performed  a  linear  least-squares  fit  of 
the  channel  number  and  energy  data  to  determine  the  pulse- 
height  axis  calibration. 

Next  we  subtracted  the  background  distribution  from  the 
data  distribution  channel  by  channel.  However,  since  the 
pulse-height  calibrations  for  corresponding  data  and 
background  files  were  usually  slightly  different,  we  first 
rebinned  the  background  distribution  into  the  appropriate 
data  distribution  bin  values  and  widths.  This  process  was 
visually  checked  for  consistency  using  peaks  in  the  data  and 
background  distributions  at  1460  and  2615  keV  due  to  the 
radioactive  nuclides  and  208ip2^^  These  nuclides  are 
found  in  ordinary  construction  materials  such  as  those  in 
the  walls  of  our  laboratory. 
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For  the  measurements  using  the  ^^Co,  ®^Sr, 

^^^Cs,  ^'^Mn,  and  sources,  the  only  further  analysis 

required  to  determine  the  response  functions  was  photopeak 
area  determination  and  subsequent  normalization  of  each 
pulse-height  distribution  to  one  incident  photon.  However, 
the  remaining  sources  emit  gamma  rays  corresponding  to  more 
than  one  energy,  or  secondary  reactions  produce  various 
energy  gamma  rays.  (See  Table  1.)  Actually,  the  ^^Co  and 
1^0*  sources  produce  multiple-energy  gamma  rays.  But  the 
two  photopeaks  at  122  and  136.5  keV  in  the  ^^Co  pulse-height 
distribution  are  not  resolved  because  of  the  BGO  resolution 
and  are  treated  as  a  single-energy  gamma  ray  corresponding 
to  123.6  keV.  This  energy  was  computed  using  a  weighted 
average  based  on  the  relative  intensities  for  the  122  and 
136.5  keV  gamma  rays.  The  6920-  and  7120-keV  gamma  rays 
emitted  from  can  be  neglected  because  of  their 

extremely  low  intensities — less  than  3  percent  of  the  6140- 
keV  gamma-ray  intensity  based  on  our  measurements.  The  511- 
keV  contribution  is  discussed  later. 

For  the  multiple-energy  sources,  the  pulse-height 
distributions  due  to  the  low-energy  gamma  rays  must  be 
subtracted  from  the  overall  multiple  gamma-ray  distribution 
in  order  to  obtain  a  response  function.  The  procedure  was 
as  follows.  We  used  the  ^^Zn  source  to  obtain  the  pulse- 
height  distribution  due  only  to  1115.5-keV  gamma  rays  by 
subtracting  the  distribution  resulting  from  the  511-keV 
gamma-rays.  The  subtraction  process  used  the  ®^Sr 
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distribution,  which  results  from  monoenergetic  514-keV  gamma 
rays.  We  multiplied  the  slope  and  intercept  of  the  ®^Sr 
pulse-height  axis  by  the  ratio  511/514  so  that  the  514-keV 
distribution  corresponded  to  the  511-keV  distribution  along 
the  pulse-height  axis.  The  appropriate  scale  factor  to 
match  the  number  of  counts  in  the  511-  and  514-keV 
photopeaks  was  determined  iteratively  by  viewing  the  Compton 
region  of  the  resultant  1115.5-keV  pulse-height 
distribution.  This  Compton  region  is  flat  when  the  511-keV 
distribution  is  scaled  appropriately.  This  same  procedure 
was  followed  for  the  ^^Na  and  multigamma  distributions. 
However,  for  the  834.8-keV  spectrum  from  ^^Mn  was  used 

in  the  subtraction  process. 

For  the  ^'^Na  source,  the  1368.6-keV  piece  of  the  pulse- 
height  distribution  was  eliminated  using  the  resultant 
1274.5-keV  spectrum  from  ^^Na.  However,  a  simple  scaling  of 
the  slope  and  intercept  of  the  1274.5-keV  spectrum  to 
approximate  the  1368.6-keV  distribution  could  not  be  used. 
Since  the  Compton  region  of  a  pulse-height  distribution  is 
relatively  flat  and  structureless  for  low-energy  gamma  rays, 
we  stretched  the  middle  of  this  region  for  1274.5  keV  so 
that  the  1274.5-keV  photopeak  was  centered  on  the  channel 
number  corresponding  to  1368.6  keV.  Furthermore,  this 
technique  preserved  the  position  of  the  backscatter  peak 
which  varies  little  with  gamma-ray  energy  for  energies 
greater  than  approximately  1000  keV.  An  analogous  process 
was  followed  to  eliminate  the  2223-keV  piece  of  the  2^®Pu/Be 
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pulse-height  distribution  using  the  resultant  2754-keV 
spectrum  and  contracting  the  Compton  region  instead  of 
stretching  it. 

The  only  other  significant  extraneous  pieces  in  the 
238pu/Be  distribution  to  be  removed  were  contributions  from 
511-  and  842-keV  gamma  rays.  We  determined  this  from 
looking  at  the  ^^^Pu/Be  pulse-height  distribution  and  also 
comparing  with  the  distribution  in  the  low  pulse-height 

regions.  The  842-keV  contribution  was  easily  removed  using 
the  ^^Mn  distribution  in  the  manner  described  previously. 
However,  the  511-keV  piece  could  not  be  removed  because  it 
is  the  result  of  annihilation  radiation  from  both  the 
extended  source  and  the  detector  assembly  materials  and 
materials  of  the  water  tank  and  shielding  housing  the 
238pu/Be  source.  Likewise,  the  511-keV  contributions  to  the 

and  ^2c*  puise-height  distributions  could  not  be 
eliminated.  Therefore,  the  response  functions  for  gamma-ray 
energies  of  4440  keV,  6140  keV,  and  11.67  MeV  should  not  be 
used  for  deconvolution  of  measured  gamma-ray  spectra  in 
pulse-height  regions  below  approximately  600  keV.  This  is 
not  a  significant  limitation  in  their  use,  however,  since 
most  gamma-ray  spectra  are  complicated  in  this  region  by 
annihilation  radiation  and  useful  information  is  difficult 
to  extract. 

Lastly,  for  those  pulse-height  distributions  with  low 
statistics  and  fairly  large  contributions  from  extraneous 
gamma  rays,  we  applied  a  five-point  data  smoothing  procedure 
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to  that  part  of  the  distribution  affected  by  the  subtraction 
process.  Figure  5  shows  the  pulse-height  distributions 
taken  with  the  2.54-cin  x  2.54-cin  BGO  detector  from  the  ®^Zn, 
SSy^  ^^Na,  and  ^^^Pu/Be  multiple  gamma-ray  sources  and  the 
distributions  after  subtraction  of  the  lower-energy  gamma- 
ray  contributions  and  subsequent  smoothing.  Figure  6  shows 
the  equivalent  pulse-height  distributions  for  the  7.62-cm  x 
7.62-cm  BGO  detector. 

A  further  correction  was  applied  to  the  2754-keV 
spectrum,  resulting  from  ^^Na,  to  correct  for  shape 
variations  with  source-to-detector  distance.  Figure  7  shows 
the  unnormalized  pulse-height  distributions  from 
acquired  with  the  2.54-cm  x  2.54-cm  detector  at  source-to- 
detector  distances  of  5,  25,  and  100  cm.  The  greatest  shape 
differences  occur  between  the  5-  and  100-cm  spectra  in  the 
pulse-height  region  below  approximately  700  keV.  This  is 
also  true  for  the  spectra  acquired  with  the  7.62-cm  x 

7.62-cm  BGO  detector.  Therefore,  the  2754-keV  distribution 
was  modified  below  700  keV  using  the  1836-keV  distribution 
taken  at  100  cm.  Due  to  the  small  differences  between  the 
®®Y  spectra  for  25-  and  100-cm  source-to-detector  distances, 
we  did  not  correct  the  and  distributions  since 

they  were  taken  at  source-to-detector  distances  of  34  cm. 

The  final  step  in  obtaining  the  response  functions  was 
to  normalize  the  single-energy  pulse-height  distributions  to 
one  gamma  ray  incident  on  the  detector.  For  the  gamma-ray 
sources  of  known  intensities,  this  was  straightforward. 
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Figure  5.  Pulse-height  distributions  acquired  with  a  2.54- 
cm  X  2.54-cin  BGO  detector  for  gamma  rays  from  ®^Zn,  °°Y , 
^^Na,  and  ^^®Pu/Be  radioactive  sources.  The  dashed  curves 
include  the  effects  of  all  gamma-rays  emitted  from  the 
sources.  The  solid  curves  are  the  distributions  from  only 
the  highest  energy  gamma  rays  emitted  from  the  sources.  The 
pulse-height  distributions  are  not  normalized. 
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Figure  6.  Pulse-height  distributions  acquired  with  a  7.62- 
cm  X  7.62-cin  BGO  detector  for  gamma  rays  from  ^^Zn,  ®®Y, 
^^Na,  and  ^^®Pu/Be  radioactive  sources.  The  dashed  curves 
include  the  effects  of  all  gamma-rays  emitted  from  the 
sources.  The  solid  curves  are  the  distributions  from  only 
the  highest  energy  gamma  rays  emitted  from  the  sources.  The 
pulse-height  distributions  are  not  normalized. 
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Pulse  Height  (keV) 


Figure  7.  Pulse-height  distributions  acquired  with  a  2.54- 
cm  X  2.54-citi  BGO  detector  for  gamma  rays  from  a 
radioactive  source  at  various  source-to-detector  distances: 

5  cm  (dashed  curve),  25  cm  (dotted  curve),  and  100  cm  (solid 
curve).  The  distributions  are  normalized  relative  to  each 
other  based  on  the  1836-keV  photopeak. 


involving  only  the  calculation  of  the  fractional  solid  angle 
subtended  by  a  detector  at  the  source  and  correcting  for 
source  decay.  For  the  remaining  sources,  we  normalized  the 
single-energy  pulse-height  distributions  for  the  7.62-cm  x 
7.62-cm  detector  using  the  average  intrinsic  photopeak 
efficiencies  for  a  7.62-cm  x  7.62-cm  BGO  detector  with  a 
source-to-detector  distance  of  91.7  cm  reported  in  Ref.  17. 
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Table  4  compares  the  intrinsic  photopeak  efficiencies  of 
Ref.  17  to  those  measured  for  our  7.62-cm  x  7.62-cm 
detector.  The  photopeak  areas  needed  to  obtain  the 
photopeak  efficiencies  in  Table  4  were  determined  by  fitting 
Gaussian  functions  to  the  data  points  above  the  half-maximum 
peak  height  values  using  MINUIT  (Ref.  20).  The  agreement  is 
good,  especially  at  the  higher  gamma-ray  energies  and 
justifies  the  use  of  Ref.  17  for  normalizing  the 
distributions  from  the  sources  of  unknown  strength.  The 
errors  shown  in  the  table  result  from  a  maximum  4.6  percent 
error  in  the  known  source  intensities  and  a  2  percent  error 
in  photopeak  areas. 

Normalization  of  the  pulse-height  distributions  from 
the  uncalibrated  sources  for  the  2.54-cm  x  2.54-cm  BGO 
detector  using  the  efficiencies  of  Ref.  17  was  accomplished 
as  follows.  We  determined  the  photopeak  areas  of  the 
uncalibrated  source  distributions  for  the  7.62-cm  x  7.62-cm 
detector.  Also,  since  these  sources  were  extended  sources, 
we  performed  Monte  Carlo  calculations  to  determine  the 
various  fractional  solid  angles  subtended  by  the  7.62-cm  x 
7.62-cm  detector  at  the  sources.  Using  the  fitted  photopeak 
areas,  computed  solid  angles,  and  the  intrinsic  efficiencies 
of  Ref.  17,  we  could  then  determine  the  strength  of  the 
2^Na,  2^®Pu/Be,  ^®0*,  and  sources.  We  then  used  these 

strengths  along  with  computed  fractional  solid  angles  and 
fitted  photopeak  areas  pertaining  to  the  2.54-cm  x  2.54-cm 
detector  to  normalize  the  single-energy  pulse-height 
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Table  4.  Instrinsic  Photopeak  Efficiencies  At  Various 


Gamma-Ray  Energies 
Detector. 

For  A  7.62-cm  X  7 

.62-cm  BGO 

Scintillation 

Gamma-Ray  Energy 
(keV) 

Instrinsic 
Ref.  17 

Photopeak 

This 

Efficiency 

Work 

123.6 

0.75 

0.89 

±0.05 

279.2 

0.95 

0.82 

±  0.04 

514 

0.90 

0.85 

±  0.03 

661.6 

0.84 

0.77 

±  0.03 

834.8 

0.79 

0.81 

±  0.04 

1115.5 

0.71 

0.71 

±  0.03 

1274.5 

0.67 

0.64 

±  0.03 

1836 

0.57 

0.61 

±  0.03 

Source-to-detector  distance  of  91.7  cm. 
Source-to-detector  distance  of  100  cm. 


distributions  from  these  sources  for  the  2.54-cm  x  2.54-cm 
detector.  The  resultant  intrinsic  photopeak  efficiencies 
for  this  detector  are  shown  in  Table  5.  For  the  lowest 
gamma-ray  energies  through  1836  keV,  the  errors  are  a 
combination  of  a  maximum  4.6  percent  error  in  the  known 
source  intensities  and  a  2  percent  error  in  photopeak  areas. 
The  remaining  efficiencies  include  a  2  percent  error  in 
fractional  solid  angle,  a  2  percent  error  in  photopeak 
areas,  and  an  estimated  5  percent  error  due  to  the 
normalization  procedure.  No  intrinsic  photopeak  efficiency 
could  be  determined  for  the  gamma-ray  energy  of  11.67  MeV 
because  a  photopeak  was  not  observed  due  to  the  low 
efficiency  at  this  energy  and  the  resolution  of  the 
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Table  5.  Instrinsic  Photopeak  Efficiencies  At  Various 
Gamma-Ray  Energies  For  A  2.54-cm  x  2.54-cm  BGO  Scintillation 
Detector  With  A  Source-To-Detector  Distance  Of  100  cm. 


Gamma-Ray  Energy 
(keV) 

Instrinsic  Photopeak  Efficiency 

123.6 

0.88 

+ 

0.05 

279.2 

0.82 

+ 

0.04 

514 

0.67 

+ 

0.02 

661.6 

0.52 

+ 

0.02 

834.8 

0.47 

+ 

0.02 

1115.5 

0.33 

+ 

0.01 

1274.5 

0.28 

+ 

0.01 

1836 

0.22 

+ 

0.01 

2754 

0.15 

+ 

0.01 

4440 

0.078 

:  0.004 

6140 

0.060 

;  0.003 

detector.  However,  the  ll.67-MeV  pulse-height  distribution 
is  correctly  normalized  based  on  the  photopeak  area  of  the 
4440-keV  peak  since  the  source  emits  4440-keV  and 

11.67-MeV  gamma  rays  of  equal  intensity. 
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IV.  RESULTS  AND  DISCUSSION 


Detector  Energy  Resolutions 

An  important  parameter  for  the  use  of  any  scintillation 
detector  is  the  percentage  energy  resolution  as  a  function 
of  gamma-ray  energy.  Figures  8  and  9,  respectively,  show 
the  2.54-cm  X  2.54-cm  and  7.62-cm  x  7.62-cm  BGO  detector 
percentage  resolutions  as  a  function  of  gamma-ray  energy. 

The  resolutions  were  calculated  for  the  full  widths  at  half 
maximums  of  Gaussian-f unction  fits  to  the  photopeaks.  The 
solid  lines  are  least-squares  fits  to  the  data  points  and 
are  given  by 

2.54-cm  X  2.54-cm  BGO  detector: 

resolution  (percent)  =  221. 524E'-y(keV)”0. 45135 
7.62-cm  X  7.62-cm  BGO  detector: 

resolution  (percent)  =  221.863E'y(keV)“0*^2863^ 

A  few  of  the  measured  resolution  values  for  the  2.54-cm  x 
2.54-cm  detector  are  11.6  ±  0.2  percent,  7.2  ±  0.1  percent, 
and  4.5  ±  0.1  percent  for  gamma-ray  energies  of  661.6,  1836, 
and  6140  keV,  respectively.  The  errors  are  2  percent  errors 
estimated  from  the  variation  in  fits  to  the  photopeaks. 

These  resolution  values  agree  well  with  the  values  reported 
in  Ref.  13:  11.8  percent,  7.2  percent,  and  4  percent.  The 

corresponding  measured  resolution  values  for  the  7.62-cm  x 
7.62-cm  detector  are  13.0  ±  0.3  percent,  8.5  ±  0.2  percent, 
and  5.6  ±  0.1  percent.  Again,  these  values  agree  well  with 
other  reported  values  (Refs.  8,10,15,17),  and  in  many  cases 
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Figure  8.  The  percentage  energy  resolution  (FWHM)  as  a 
function  of  gamina-ray  energy  for  the  2.54-cin  x  2.54-cin  BGO 
scintillation  detector. 
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Figure  9.  The  percentage  energy  resolution  (FWHM)  as  a 
function  of  gamma-ray  energy  for  the  7.62-cm  x  7.62-cm  BGO 
scintillation  detector. 


our  7.62-cm  X  7.62-cm  BGO  detector  exhibits  slightly  better 
resolution. 

Response  Functions 

The  response  function  curves,  plots  of  the  number  of 
counts  recorded  by  the  detectors  per  unit  pulse  height  per 
incident  gamma  ray,  are  shown  in  Figs.  10  and  11  for  the 
2.54-cm  X  2.54-cm  BGO  detector  and  Figs.  12  and  13  for  the 
7.62-cm  X  7.62-cm  BGO  detector.  Some  general  features  are 
the  following.  The  backscatter  peak  occurring  at  a  pulse 
height  of  approximately  200  to  250  keV  is  a  prominent 
feature  of  every  response  function  for  both  detectors.  The 
size  of  this  peak  in  relation  to  the  photopeak  generally 
increases  as  the  energy  of  the  photopeak  increases.  Since 
the  detectors  were  well  removed  from  the  walls,  the  ceiling, 
the  floor,  and  any  apparatus  in  the  laboratory  room  during 
most  of  the  response  function  measurements,  the  backscatter 
peaks  must  be  primarily  the  result  of  Compton-scattered 
gamma  rays  from  the  detector  housing  materials  and  the 
aluminum  camera  tripods.  Thus,  the  relative  importance  of 
the  backscatter  peak  to  a  given  response  function  will 
depend  on  the  construction  of  the  detector,  its  mounting, 
and  other  nearby  materials  and  will  obviously  vary  widely 
between  detectors  and  experiments.  A  similar  concern 
applies  to  the  annihilation  radiation  occurring  at  a  pulse 
height  of  511  keV  in  all  the  response  functions  for  E-y  > 
1836  keV.  The  annihilation  peak  arises  from  gamma-ray  pair 
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Figure  10.  Response  function  curves  for  the  2.54-cm  x  2.54 
cm  BGO  scintillation  detector  at  gamma-ray  energies  of 
123.6,  279.2,  514,  661.6,  834.8,  and  1115.5  keV. 
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Figure  11.  Response  function  curves  for  the  2.54-cin  x  2.54- 
cm  BGO  scintillation  detector  at  gamma-ray  energies  of 
1274.5,  1836,  2754,  4440,  6140,  and  11,670  keV. 
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Figure  12.  Response  function  curves  for  the  7.62-cin  x  7.62 
cm  BGO  scintillation  detector  at  gam^a-ray  energies  of 
123.6,  279.2,  514,  661.6,  834.8,  and  1115.5  keV. 
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Figure  13.  Response  function  curves  for  the  7.62-cin  x  7.62- 
cm  BGO  scintillation  detector  at  gamma-ray  energies  of 
1274.5,  1836,  2754,  4440,  6140,  and  11,670  keV. 
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production  within  materials  surrounding  the  detector  and 
again  is  specific  to  the  detector  and  experiment. 

Therefore,  the  pulse-height  regions  below  approximately  600 
keV  of  the  response  functions  shown  in  Figs.  10-13  are 
detector  specific.  When  using  these  response  functions  to 
unfold  gamma-ray  spectra,  the  unfolded  results  below  600  keV 
should  be  used  with  caution. 

For  the  2.54-cm  x  2.54-cm  BGO  detector,  the  single-  and 
double-escape  peaks  become  prominent  at  E-y  =  2754  keV  and 
dominate  the  response  at  E-y  >  4440  keV.  However,  the 
photopeak  is  always  the  dominant  part  of  the  response 
function  for  the  7.62-cm  x  7.62-cm  BGO  detector.  Also,  the 
Compton  regions  of  the  7.62-cm  x  7.62-cm  responses  are 
significantly  less  important  than  these  regions  for  the 
2.54-cm  X  2.54-cm  detector.  Therefore,  if  energy  resolution 
is  not  a  primary  concern,  the  larger  BGO  detector  is  a 
better  choice  for  measurement  of  continuous  gamma-ray 
spectra  since  most  of  its  response  is  in  the  photopeak.  A 
photopeak-dominated  response  can  be  used  to  simplify  the 
unfolding  process  and,  at  the  least,  results  in  better 
unfolded  results  using  response  matrix  unfolding  methods. 
However,  the  energy  resolution  of  the  2.54-cm  x  2.54-cm 
detector  is  better  than  the  7.62-cm  x  7.62-cm  detector  and 
should  be  used  when  resolution  is  important.  Also,  its 
smaller  size  makes  it  a  better  choice  when  measuring  high- 
energy  gamma  rays  accompanied  by  large  numbers  of  neutrons 
and  low-energy  gamma  rays. 
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V.  RESPONSE  MATRIX  GENERATION 


Technique 

In  order  to  use  the  response  functions  reported  in  this 
paper  for  unfolding  measured  gamma-ray  spectra,  most 
unfolding  methods  require  a  response  matrix.  The  general 
response  matrix  is  a  square  matrix  of  elements  with  each 
element  denoted  by  Rj]^(P  j  ,E]^)  .  Rjj^(Pj,Ej^)  is  the 
probability  that  a  gamma  ray  incident  on  the  detector  with 
an  energy  in  the  bin  Ej^  of  width  produces  a  count  in 

the  pulse-height  distribution  in  the  bin  Pj  of  with  Z5iPj. 
Rjj^(Pj,E)^)  is  related  to  the  measured  gamma-ray  spectrum  and 
incident  gamma-ray  spectrum  through  the  relation 
Mj(Pj)  =  p3„(Pj,EK)T„(Ek), 

where  the  T)^(E]^)'s  compose  a  discrete  representation  of  the 
gamma-ray  spectrum  incident  on  the  detector  and  the  Mj(Pj)'s 
are  a  discrete  representation  of  the  measured  spectrum.  To 
determine  the  incident  gamma-ray  spectrum,  T(E),  the 
response  elements,  R  jj^  ( P  j  ,  Ej^ )  ,  are  needed  at  many  energy  and 
pulse-height  values  for  arbitrary  Z!iPj  and  bin  widths. 

We  measured  the  response  functions  for  2.54-cm  x  2.54-cm  and 
7.62-cm  X  7.62-cm  BGO  detectors  at  only  12  energy  values 
spanning  the  gamma-ray  energy  range  of  123.6  keV  to  11.67 
MeV.  However,  using  an  interpolation  procedure  discussed 
below,  these  12  responses  can  be  used  to  obtain  response 
functions  at  any  values  of  gamma-ray  energy  between  123.6 
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keV  and  11.67  MeV  and  thus  for  construction  of  a  response 
matrix. 

Generation  of  a  response  function  for  a  gamma-ray 
energy  not  measured  is  done  as  follows.  Suppose  that  one 
column  of  the  response  matrix  is  the  energy  bin  centered  at 
Ej^  =  5290  keV.  This  energy  is  bracketed  by  the  energies  of 
4440  keV  and  6140  keV  for  which  measured  response  functions 
exist.  The  5290-keV  response  function  is  interpolated  from 
the  4440-keV  and  6140-keV  response  functions.  First,  the 
4440-keV  response  function  curve  is  shifted  and  stretched 
along  the  pulse-height  axis  so  that  the  positions  of  the 
photopeak  (4440  keV) ,  single-escape  peak  (3929  keV) ,  and 
double-escape  peak  (3418)  acquire  new  positions 
corresponding  to  the  positions  of  these  three  peaks  for  a 
5290-keV  response  function — 5290  keV,  4779  keV,  and  4268 
keV,  respectively.  The  magnitude  of  the  4440-keV  response 
curve  is  not  modified.  Also,  the  positions  of  the 
backscatter  and  annihilation  peaks  in  the  shifted  and 
stretched  4440-keV  response  curve  are  left  unchanged.  The 
above  is  accomplished  by  dividing  the  pulse-height  axis  of 
the  4440-keV  response  curve  into  three  regions:  pulse 
heights  <  733  keV,  pulse  heights  >  3167  keV,  and  733  <  pulse 
heights  <  3167  keV.  These  regions  are  for  the  2.54-cm  x 
2.54-cm  BGO  response  function  and  are  different  for  the 
7.62-cm  X  7.62-cm  detector  because  of  the  difference  in 
energy  resolution.  The  transformation  from  the  4440-keV 
response  to  the  shifted  and  stretched  4440-keV  response 
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consists  of  leaving  the  pulse  heights  <  733  keV  untouched, 
adding  a  constant  value  of  850  keV  (5290  -  4440)  to  all 
pulse  heights  >  3167  keV,  and  multiplying  the  middle  region 
of  pulse  heights  by  the  appropriate  stretching  factor  to 
fill  the  middle  region  of  the  shifted  and  stretched  4440-keV 
response  pulse-hej.ght  axis.  This  procedure  maintains  the 
shapes  of  all  peaks  in  the  response  function.  Second,  the 
6140-keV  response  function  curve  is  shifted  and  compressed 
in  the  manner  described  above  except  that  the  shift  is 
negative  and  the  three  pulse-height  regions  are  changed  so 
as  to  correspond  to  the  shape  of  the  6140-keV  response. 

After  this  shifting  and  stretching  or  compressing  procedure, 
the  4440-keV  and  6140-keV  pulse-height  axes  appear  as  5290- 
keV  pulse-height  axes.  The  final  step  is  a  weighted 
logarithmic  interpolation  of  the  ordinate  axes  of  the 
modified  4440-keV  and  6140-keV  response  curves  on  a  point- 
by-point  basis  along  the  common  pulse-height  axis  to 
determine  the  5290-keV  curve.  The  weighting  is  based  on  the 
natural  logarithm  of  the  energy  differences  between  the 
three  response  functions.  The  interpolation  procedure  just 
described  is  shown  graphically  in  Fig.  14  for  the  2.54-cm  x 
2.54-cm  BGO  detector. 

To  test  this  construction  procedure  for  response 
functions,  response  functions  for  both  detectors  were 
generated  at  gamma-ray  energies  for  which  measured  response 
functions  existed.  Two  examples  for  the  2.54-cm  x  2.54-cm 
BGO  detector  at  E-y  =  1274.5  and  4440  keV  are  shown  in  Figs. 
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Figure  14.  (a)  Resp'^nse  function  curve  for  the  2.54-cin  x 

2.54-cin  BGO  scintillation  detector  at  E-r  =  4440  keV  in 
which  the  pulse-height  axis  has  been  shifted  and  stretched 
to  match  the  pulse-height  axis  for  E-y  =  5290  keV.  (b) 
Response  function  curve  at  E-y  =  6140  keV  in  which  the 
pulse-height  axis  has  been  shifted  and  compressed  to  match 
the  pulse-height  axis  for  E-y  =  5290  keV.  (c)  Response 
function  curve  at  E-y  =  5290  keV  interpolated  from  the  4440- 
and  6140-keV  response  curves. 
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15  and  16.  The  interpolated  1274.5-keV  response  used  the 
measured  1115.5-keV  and  1836-keV  response  functions;  the 
interpolated  4440-keV  response  used  the  measured  2754-keV 
and  6140-keV  response  functions.  There  is  very  good 
agreement  between  the  interpolated  and  measured  response 
curves,  especially  for  pulse  heights  greater  than  600  keV 
where  the  responses  are  applicable. 

The  interpolated  response  functions  are  the  basis  for 
obtaining  the  response  matrix.  A  column  of  the  matrix 
Rjj^(Pj,Ej^)  is  the  appropriately  binned  (in  pulse  height) 
response  function  for  a  gamma-ray  energy  of  and  energy 
bin  width  of  We  chose  to  construct  this  response 

function  as  a  simple  average  of  three  interpolated  response 
functions  at  energy  values  equally  spaced  within  A.Ej^ 
centered  at  Ej^.  The  averaging  process  accounts  for  possible 
variations  in  the  detector  response  function  across  the 
energy  bin  AEj^.  These  variations  are  most  important  if 
2^Ej^  is  large.  Also,  the  averaging  process  effectively 
smooths  the  response  function,  which  is  beneficial  in  the 
unfolding  process. 

Computer  Code 

The  computer  code,  BGRESP,  used  to  construct  the 
response  matrix  for  a  2.54-cm  x  2.54-cm  or  7.62-cm  x  7.62-cm 
BGO  detector  using  the  experimental  response  functions 
reported  in  this  paper  is  listed  in  the  Appendix.  The  code 
is  written  in  FORTRAN-77.  The  primary  components  of  the 
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Figure  15.  Measured  (dashed  curve)  and  interpolated  (solid 
curve)  response  functions  for  the  2.54-cm  x  2.54-cni  BGO 
scintillation  detector  at  E-y  =  1274.5  keV. 


Pulse  Height  (keV) 

Figure  16.  Measured  (dashed  curve)  and  interpolated  (solid 
curve)  response  functions  for  the  2.54-cin  x  2.54-cm  BGO 
scintillation  detector  at  E-y  =  4440  keV. 


code  are  the  main  program  and  two  subroutines.  The  two 
subroutines  are  identical  except  for  the  initial  data 
statements.  One  subroutine  constructs  an  interpolated 
response  function  for  the  2.54-cm  x  2.54-cm  detector  at  a 
specific  value  of  gamma-ray  energy  passed  by  the  main 
program;  the  other  subroutine  constructs  interpolated 
response  functions  for  the  7.62-cra  x  7.62-cm  detector.  The 
main  program  forms  the  response  matrix  from  the  interpolated 
response  functions.  The  maximum  size  allowed  for  the 
response  matrix  is  100  gamma-ray  energy  bins  by  100  pulse- 
height  bins. 

The  input  and  output  of  the  code  is  simple.  The  input 
consists  of  specifying  one  of  the  two  sizes  of  BGO 
detectors,  2.54-cm  x  2.54-cm  or  7.62-cm  x  7.62-cm,  the  range 
of  pulse  heights  for  generation  of  the  response  matrix,  and 
the  energy  and  pulse-height  bin  width  of  a  response  matrix 
element.  This  bin  width  is  entered  as  a  multiple  of  the 
detector  energy  resolution.  Thus  the  response  matrix 
contains  variable  bin  widths  which  increase  as  the  energy 
increases.  The  output  c/onsists  of  three  files:  (1)  a 
listing  of  the  response  matrix  in  an  easily  readable  format, 
(2)  a  listing  of  the  response  matrix  in  a  compact  format  for 
input  to  an  unfolding  code,  and  (3)  a  listing  of  the  bin 
values  of  the  response  matrix  for  input  to  an  unfolding 
code. 

Several  tests  of  the  response  matrix  computer  code  were 
performed  for  both  BGO  detectors.  Fig.  17  is  an 
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illustration  of  one  such  test.  Fig.  17a  is  a  simulated 
pulse-height  distribution  from  the  2.54-cm  x  2.54-cm 
detector  due  to  six  equal-intensity  gamma-ray  sources 
obtained  by  adding  together  the  measured  response  functions 
at  E-y  =  834.8,  1274.5,  1836,  2754,  4440,  and  6140  keV. 

These  six  response  functions  were  scaled  before  adding  them 
together  so  that  each  one  corresponds  to  the  same  number  of 
gamma  rays  incident  on  the  detector.  Furthermore,  the 
resultant  pulse-height  distribution  was  normalized  to 
correspond  to  a  total  of  six  gamma  rays  incident  on  the 
detector.  The  fact  that  the  simulated  pulse-height 
distribution  is  a  sum  of  measured  response  functions  instead 
of  an  actual  simultaneous  measurement  of  six  gamma-ray 
sources  is  not  a  limitation  of  the  test  since  the  response 
matrix  results  from  an  averaging  of  the  interpolated 
response  functions  and  not  directly  from  the  measured 
response  functions.  Unfolding  this  pulse-height 
distribution  with  the  2.54-cm  x  2.54-cm  response  matrix 
should  ideally  yield  a  gamma-ray  spectrum,  T(E),  consisting 
of  six  delta-function  peaks  at  the  energies  834.8,  1274.5, 
1836,  2754,  4440,  and  6140  keV,  each  peak  with  an  area  of 
one.  Fig.  17b  shows  the  unfolded  spectrum  T(E).  A  response 
matrix  of  size  92  x  92  with  a  bin  wid;:h  of  0.5  times  the 
detector  resolution,  generated  from  the  code  BGRESP,  was 
used  for  the  unfolding.  The  areas  of  the  peaks  are  1.03, 
0.96,  0.89,  1.00,  1.02,  and  1.03,  for  a  total  area  of  5.93 
as  compared  to  the  ideal  expected  total  area  of  6.  A  test 
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Figure  17.  (a)  A  simulated,  measured  pulse-height 

distribution  resulting  from  equal-intensity  gamma-ray 
sources  at  E-y  =  834.8,  1274.5,  1836,  2754,  4440,  and  6140 
keV  for  the  2.54-cm  x  2.54-cm  BGO  scintillation  detector. 

(b)  The  unfolded  gamma-ray  spectrum  resulting  from  (a)  using 
a  92  X  92  response  matrix  constructed  from  measured  2.54-cm 
X  2.54-cm  response  functions. 
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of  the  response  matrix  computer  code  for  the  7. 6 2 -cm  x  7.62- 
cm  detector  is  illustra+'ed  in  Fig.  1 8 .  Again,  Fig.  18a 

is  a  simulated  pulse -height  distribution  from  the  7.62-cm  x 
7.62-cm  detector  due  to  six  equal-intensity  gamma-ray 
sources.  Fig.  18b  shows  the  unfolded  spectrum.  A  response 
matrix  of  size  77  x  77  with  a  bin  width  of  0.5  times  the 
detector  resolution  was  used  for  the  unfolding.  The  areas 
of  the  peaks  are  0.996,  1.02,  1.04,  1.04,  l.oo,  1.04,  for  a 
total  area  of  6.14  as  compared  to  the  ideal  expected  total 
area  of  6. 
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Figure  18.  (a)  A  simulated,  measured  pulse-height 

distribution  resulting  from  equal-intensity  gamma-ray 
sources  at  E-y  =  834.8,  1274.5,  1836,  2754,  4440,  and  6140 
keV  for  the  7.62-cm  x  7.62-cm  BGO  scintillation  detector. 

(b)  The  unfolded  gamma-ray  spectrum  resulting  from  (a)  using 
a  77  X  77  response  matrix  constructed  from  measured  7.62-cm 
X  7.62-cm  response  functions. 
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VI .  CONCLUSIONS 


The  24  response  functions  spanning  the  gaituna-ray  energy 
range  of  123.6  keV  to  11.67  MeV  for  2.54-CTn  x  2.54-cin  and 
7.62-cin  X  7.62-cin  BGO  scintillation  detectors  along  with  the 
computer  code  for  calculating  response  matrices  should  be 
useful  by  others  possessing  the  same  size  detectors  with 
equivalent  or  worse  detector  resolutions.  They  can  be  used 
for  unfolding  gamma-ray  pulse-height  distributions  with 
reliable  results  for  pulse  heights  above  600  keV.  Since  the 
response  functions  correspond  to  a  source-to-detector 
distance  of  100  cm,  they  can  be  used  with  little  error  for 
any  experimental  arrangement  in  which  the  BGO  detector  is  at 
least  approximately  50  cm  from  the  source. 
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APPENDIX 


The  following  computer  code,  BGRESP,  calculates  a 
response  matrix  for  a  2.54-cm  x  2.54-cm  or  7.62-cm  x  7.62-cm 
BGO  detector  using  a  set  of  12  experimental  response 
functions  for  each  detector  covering  the  gamma-ray  energy 
range  of  123.6  keV  to  11.67  MeV.  A  copy  of  this  code  along 
with  the  24  response  function  data  files  and  sample  input 
and  output  can  be  obtained  from  the  authors  by  writing  to 
their  address  or  calling  the  number  shown  in  the  listing  of 
the  code. 


PROGRAM  BGRESP 
C 

C  by: 

C  Rex  R.  Kiziah 

C  Physics  Department 

C  United  States  Air  Force  Academy 

C  Colorado  Springs,  CO  80840-5701 

C  (719)  472-3355 

C 

C  VERSION  1.0,  3  April  1990 
C 

C  Program  BGRESP  is  a  FORTRAN  code  to  calculate  a  gamma-ray  response 
C  matrix  for  a  2.54-cm  x  2.54-cm  or  a  7.62-cn  x  7.62-cm  BGO  inorganic 
C  scintillator.  The  maximum  matrix  size  is  100  x  100  and  can  span  any 
C  enerjy  and  puls-'-'ieight  range  between  123.6  and  11670.0  lieV.  The 
C  matrix  is  in  a  format  such  that  the  columns  are  pulse-height  bins  and 
C  the  rows  are  energy  bins  (each  in  l?eV).  The  2-D  array  FIHRESP(I,J) 

C  contains  the  matrix  (1=1, 100, J-l, 100  max).  F1NRESP(1',J')  is  the 
C  probability  for  one  gamma  ray  in  the  energy  bin  1'  to  produce  a  count 
C  in  the  pulse-height  bin  J'.  The  matrix  is  calculated  via  interpolation 
C  from  12  data  files  for  each  detector.  Each  data  file  is  an  experimental 
C  response  function— bacliground  corrected,  pulse  height  calibrated, 

C  properly  normalized,  etc.  The  response  functions  were  acquired  for  a 
C  2.54-cm  X  2.54-cn  and  a  7.62-cm  x  7.62-cm  BGC  scintillator,  each 
C  purchased  from  BICROM  Corporation.  The  resolution  for  the  2.54-cm  x 
C  2.54-cm  detector  was  determined  experimentally  to  be  R  =  DeltaE/E  = 

C  2. 21524*E(gamma) **-0.45135,  with  E(gamma)  in  keV.  For  the  7.62-cm  x 
C  7.62-cm  detector,  R  =  2.21863*E(gamma)**-0.42863.  If  this  code  is  to 
C  be  used  for  a  BGO  detector  with  poorer  resolution,  the  user  should  first 
C  fold  the  data  files  with  a  Gaussian  of  the  proper  width  to  simulate 
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the  poorer  resolution  effects  and  then  appropriately  nodify  the 

resolution  fornulas  used  in  this  code. 

Running  the  Code: 

(1)  The  user  must  have  the  24  response  function  data  files,  12  for 
the  2.64-cb  X  2.54-cb  BGO  detector  and  12  for  the  7.62-cb  x 
7.62-cb  detector. 

(2)  If  the  user's  detector  resolution  is  worse  than  that  given  above 
in  the  description  of  the  code,  appropriate  Bodifications  to  the 
data  files  and  this  code  Bust  be  Bade. 

(3)  Run  BGRESP.  The  user  is  proBpted  for  all  necessary  inforaation. 

(4)  Three  output  files  are  generated.  BGRESP.OUT-'a  file  listing  the 
response  Batrix  bins  and  the  response  satrix  in  a  foraat  for  ease 
of  read:.'.;.  BGO. HAT— a  file  containing  the  response  natrix 

in  a  foraat  for  input  to  an  unfolding  code.  BGLIMS.DAT—  a  file 
containing  the  number  of  bins  and  the  response  matrix  bin 
values  (in  HeV)  for  input  to  an  unfolding  code. 

Program  Summary : 

Subroutines  Used: 

BGIRESP  ==  Calculates  a  gamma-ray  response  function  for  any 
gamma-ray  energy  in  the  range  123.6  to  11670.0  keV 
for  the  2.54-cb  x  2.54'CB  BGO  detector.  The  1-D  array 
RESP(I),  1=1000,  contains  the  response  function. 
RESP(r)  is  the  probability  per  unit  pulse  height  (in 
keV)  that  one  gamma  ray  of  the  specified  energy 
produces  a  count  in  the  pulse  height  bin  I'.  That  is, 
RESP(I')  has  units  of  Counts/keV/gamma-ray. 

BG3RESP  ==  Identical  to  BGIRESP  except  that  the  calculation  is  for 
the  7,62-cn  x  7.62-cb  BGO  detector. 

Functions  Used: 

TERP  ==  An  interpolation  routine. 

Common  Blocks  Used: 

DATARAVS  ==  Contains  2  1-D  arrays  shared  between  the  main  program 
and  the  two  subroutines. 

PARK  ==  Contains  3  parameters  shared  between  the  main  program 
BGRESP  and  the  subroutines  BGIRESP  and  BG3RESP. 

Local  Variables  Used: 

BINC  ==  Real  variable.  A  variable  for  storing  the  centroid 
of  each  bin  of  the  response  natrix.  In  keV. 

BINFAC  ==  Real  variable.  Factor  for  determining  individual 
bin  sizes  for  response  natrix.  It  is  a  multiple  of 
the  detector  energy  resolution— thus,  bin  sizes  are 
unequal. 

BINW  ==  Real  variable.  A  variable  for  storing  the  width 
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of  each  bin  of  the  response  natrix.  In  keV. 

DELTA  ==  Real  variable.  Used  in  computation  of  energy  values 
for  the  response  functions  used  in  construction  of 
the  response  matrix. 

ERESP  ==  Real  1-D  300  element  array.  Contains  the  energy  values 
(in  keV)  for  generation  of  response  functions  used  to 
construct  the  response  matrix. 

FINRESP  ==  Real  2-D  100  x  100  array.  Contains  the  response  matrix. 
FINRESP(r,J')  is  the  probability  that  a  gamma  ray  in 
the  energy  bin  r  produces  a  count  in  the  pulse-height 
bin  J'.  Both  bin  widths  are  in  keV. 

IDET  ==  Integer  variable.  Value  is  1  or  2.  1  is  for  2.54-cm 
X  2.54-cm  detector;  2  is  for  7.62-cn  x  7.62-cm. 
Determines  tor  which  detector  the  response  matrix  is 
computed. 

MH  ==  Integer  variable.  A  counter. 

NBINS  ==  Integer  variable.  The  i  of  bins  in  the  response  matrix. 

Cannot  exceed  100. 

NN  Integer  variable.  A  counter. 

NRESP  ==  Integer  variable.  This  is  the  number  of  response 
functions  used  to  compute  the  response  matrix. 
NRESP=MBINS*3.  3  response  functions  per  bin  are  used. 

PHBIN  ==  Real  variable.  In  Common  Block  PARK.  The  pulse-height 

bin  width  (in  keV)  for  each  response  function  computed 
in  BGIRESP  or  BG3RESP. 

PHLIHS  ==  Real  1-D  101  element  array.  Contains  the  pulse- 

height  values  (in  keVj  for  the  response  matrix  bins. 

PRANH  ==  Real  variable.  In  Common  Block  PARK.  Highest  value 

of  pulse  height  (in  keV)  for  generation  of  response 
matrix. 

PRANL  ==  Real  variable.  In  Common  Block  PARK.  Lowest  value 

or  pulse  height  (in  keVj  for  generation  of  response 
matrix. 

RBRESP  ==  Real  1-D  100  element  array.  Same  as  RESP  with 

the  bin  width  integrated  out,  but  rebinned  to  the 
response  matrix  bin  sizes. 

RESP  -=  Real  1-D  1000  element  irray.  In  Common  Block 

DATARAYS.  Array  to  s  ore  response  function  for  a 
'i  ’en  gamma-ray  energy  RESPd')  is  the 
probability  per  unit  p  Ise  height  (in  keV)  for  one 
gamma  ray  to  produce  a  count  in  the  pulse  height  bin 
r.  Also  stores  the  same  response  function  with  the 
bin  width  integrated  out—yields  probability  instead 
of  probability  per  unit  pulse  height. 

RESPPH  ==  Real  1-D  1000  element  array.  In  Common  Block 
DATARAYS.  Contains  the  pulse-height  bin  values 
(in  keV)  for  each  response  function  computed  by 
BGIRESP  or  BG3RESP. 

SUBRESP  ==  Real  2-D  300  x  100  array.  Contains  all  the  response 
functions,  RBRESP,  needed  for  computation  of  the 
response  matrix.  300— max  f  of  response  functions. 

SUHHI  ==  Real  variable.  Used  for  the  rebinning  of  response 
function  computed  by  BGIRESP  or  BG3RESP  to  bins 
of  the  response  matrix. 
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SUMLO  ==  Real  variable.  Used  for  the  rebinning  of  response 
function—see  SUHHI. 

SUHRESP  ==  Real  1-D  1000  eleient  array.  A  running  sun  of  RESP. 
SUMRESP( I ) =RESP( I ) +RESP( I-l ) +RESP( 1-2 ) + . . . .  Used 
for  rebinning  of  response  function—see  SUMHI. 

TYPE  ==  Character  variable.  Either  '2.54-cn  x  2.54-cii'  or 
'7.62-cb  X  7.62-ciii'. 


Set  up  the  coBnon  blocks  used  between  the  sain  prograi  and  the 
subroutines  BGIRESP  and  BG3RESP. 

COKHON/PARM/PRANL , PRANH , PHBIN 
COHHON/DATARAYS/RESP( 1000 ) ,RESPPH( 1000) 

Diiension  the  variables. 

CHARACTER  TYPE*17 

DIHENSION  PHLIHS(lOl) ,ERESP(300) ,SUKRESP(1000) ,RBRESP(100) , 

1  S0BRESP( 300,100) , FINRESP( 100,100) 

Establish  the  unit  for  interactive  input  of  data  needed  by  the  code. 

CALL  ASSIGN(1,'SYS$0UTPUT') 

Proipt  user  for  required  input:  Detector,  Response  Matrix  Size, 
Response  Matrix  Bin  Size 

WRITE(1,10) 

FORMAT (IX, 'Enter  the  detector  for  which  response  matrix  is' 

1  '  to  be  generated:'/lX,'{l)  for  1"  x  1",  (2)' 

2  '  for  3"  x  3":  ',$) 

READ(1,*)  IDET 

IF(IDET.NE.1.AHD.IDET.NE.2)  GO  TO  5 
WRITE(1,30) 

FORMAT (IX, 'Enter  the  range  of  pulse  heights  (in  keV)  for', 

1  '  generation  of  response  matrix — '/IX, 'Must  be  between', 

2  '  123.6  and  11670.0  keV:  ',$) 

READd,*)  PRANL, PRANH 

IF(PRANL.LT. 123. 6. OR. PRANH. GT.11670.0.0R.PRANL.GE. PRANH) 

1  GO  TO  25 
WRITE(1,50) 

FORMAT (IX, 'Enter  factor  for  determination  of  response', 

1  '  matrix  bin  size  as  a'/lX, 'multiple  (<  or  >  1)  of  the', 

2  '  detector  energy  resolution:  ',$) 

READd,*)  BINFAC 

Establish  for  which  detector  the  response  matrix  is  being  computed 
and  set  up  the  character  variable  TYPE  for  output  identification 
purposes. 

IFdDET.EQ.l)  THEN 

TYPE  =  '2.54-cb  X  2.54-cn' 

ELSE 

TYPE  =  '7.62-cb  X  7.62-cra' 
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ENDIF 


C 

C  Establish  the  pulse-height  values  (in  keV)  of  the  response  natrix 
C  bins.  The  bin  sizes  are  based  on  the  detector  resolution—they 
C  are  a  user-specified  nultiple,  BINFAC,  of  the  detector  resolution. 

C  Also,  count  the  nunber  of  bins  that  are  established  and  store 
C  this  nunber  in  NBINS. 

C 

NBINS  =  0 
PHLIHS(l)  =  PRANL 
DO  70  I  =  2,101 

NBINS  =  NBINS  +  1 
IFiIDET.EQ.l)  THEN 

PBLINS(I)  =  BINFAC*2.21524*PHLIHS{I-1)**0. 54865 
1  +  PHLIHS(I-1| 

ELSE 

PHLIHS(I)  =  BINfAC*2.21863*PHLIMS(I-l)**0. 57137 
1  +  PHLIKS{I-1) 

ENDIF 

IF(PHLIHS(I).GT.PRANH)  THEN 
PHLIHS(I)  =  PRAiffl 
GO  TO  75 

ENDIF 

70  CONTINUE 
75  CONTINUE 
C 

C  Check  to  Bake  sure  that  the  input  provided  by  the  user  for  the 
C  pulse-height  range  of  the  response  natrix  and  the  bin  size  are 
C  conpatible—that  is,  the  nax  size  of  the  response  natrix  is 
C  100  X  100  and  too  snail  of  a  bin  size  will  result  in  a  response 
C  natrix  of  size  100  x  100,  yet  it  will  not  span  the  entire  range 
C  of  pulse  heights. 

C 

IF{PHLIHS(NBINS+1).LT.PRANH)  THEN 
WRITE(1,80) 

80  FORNAT( IX, 'Factor  for  deternination  of  response  natrix', 

1  '  resolution  is  too  snail! '/IX, 'RESPONSE  NATRIX  SIZE', 

2  '  EXCEEDS  THE  ALLOWED  SIZE') 

GO  TO  4 

ENDIF 

C 

C  The  energy  and  pulse-height  bins  of  the  response  natrix  have  now 
C  been  established.  For  each  energy  bin,  a  pulse-height  integrated 
C  response  function  is  needed.  Since  the  energy  bins  nay  be  wide 
C  in  sone  cases,  to  use  as  nuch  of  the  experinental  data  as  possible, 

C  the  code  uses  three  pulse-height  integrated  response  functions  for 
C  each  energy  bin.  The  procedure  is  as  follows:  Divide  each  energy 
C  bin  into  6  egual  sub-bins.  Then  compute  a  pulse-height  integrated 
C  response  function  at  the  energy  value  of  the  first,  third,  and  fifth 
C  sub-bins.  Average  these  three  response  functions  to  obtain  the 
C  response  function  for  one  energy  bin  of  the  response  natrix.  The 
C  following  loop  sets  up  the  energy  values  for  conputation  of  the 
C  3*NBINS  response  functions. 


52 


NRESP  =  3*NBINS 
DO  100  I  =  2,NBIHS+1 

DELTA  =  PHLIHS(I)  -  PHLIHS(I-l) 

DELTA  =  DELTA/6.0 
DO  90  J  =  1,3 
K  =  1-2 

E1!ESP(3*K+J)  =  PHLr«S(I-l)  +  (2*J  -  1)*DELTA 
90  CONTINUE 

100  CONTINUE 
C 

C  This  loop  calls  either  the  subroutine  BGIRESP  or  BG3RESP  to  compute 
C  the  3*NBINS  pulse-heiqht  integrated  response  functions.  The  energy 
C  value  is  the  only  variable  passed  directly  to  the  subroutines.  The 
C  other  variables  needed  are  passed  through  the  connon  blocks  PARN  and 
C  DATARAYS.  Each  subroutine  returns  a  response  function  in  the  array 
C  RESP.  RESP  is  an  absolute  differential  efficiency  curve  for  the 
C  appropriate  detector— that  is,  a  curve  of  counts  per  unit  pulse  height 
C  (in  keV)  per  incident  gasma  ray.  RESP  is  returned  via  the  conmon 
C  block  DATARAYS.  Also  returned  is  the  array  RESPPH  which  contains 
C  the  1000  pulse  height  values  at  which  the  response  function  is  computed. 
C  Lastly,  the  constant  bin  width  of  the  response  function,  PHBIN  is 
C  returned  via  the  common  block  PARK.  This  bin  width  is  needed  so  that 
C  the  absolute  differential  efficiency  curve  can  be  converted  to  an 
C  absolute  efficiency  curve,  i.e.,  a  pulse-height  integrated  response 
C  function. 

C 

DO  150  I  =  1, NRESP 
IF(IDET.EQ.l)  THEN 

CALL  BG1RESP(ERESP(I)) 

ELSE 

CALL  BG3RESP(ERESP(I)) 

ENDIF 

C 

C  The  response  function  has  been  ronputed.  Now,  integrate  over  pulse 
C  height  and  store  the  pulse-height  integrated  response  back  into  the 
C  array  RESP.  RESP(I')  is  now  the  probability  that  a  gamma  ray  of 
C  the  specified  energy  for  which  RESP  was  computed  produces  a  count 
C  in  the  pulse  height  bin  I'. 

C 

DO  no  J  =  1,1000 

RESP(J)  =  RESP(J)*PHBIN 
no  CONTINUE 

C 

C  The  next  step  is  to  rebin  each  RESP  into  the  pulse  height  bins  of 
C  the  response  matrix.  RESP  is  a  pulse-height  integrated  response 
C  function  evaluated  at  1000  pulse-height  bins.  Rebin  into  NBINS 
C  bins  and  store  the  result  in  RBRESP. 

C 

SUMRESP(l)  =  RESP(l; 

DO  120  K  =  2,1000 

SUKRESP(K)  =  SUHRESP(K-l)  ^  RE3P(K) 

120  CONTINUE 

SUHLO  =  TERP(PHLIHS(li,RESPPH,SUHRESP) 

DO  130  J  =  2,NBINS+1 
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SUHHI  =  TERP(PHLIHS(J),PESPPH,SUHRESP) 

RBRESP(J-l)  =  SUHHI  -  SUHLO 
SUHLO  =  SUHHI 
130  CONTINUE 

C 

C  Load  all  NEESP  properly  binned,  pulse-heiqht  integrated  response 
C  functions  into  the  2-D  array  SUBRESP. 

C 

DO  140  K  =  1,NBINS 

SUBRESP(I,K)  =  RBRESP(K) 

140  CONTINUE 

150  CONTINUE 
J  =  0 
C 

C  There  are  three  response  functions  per  response  matrix  energy  bin. 

C  Convert  these  three  response  functions  into  one  response  function 
C  per  energy  bin  by  averaging  the  three.  Th^"  result  is  FINRESP— 

C  a  2-D  array  which  is  the  response  matrix.  FIN?.ESP(r,J')  is  the 
C  probability  that  a  gamma  ray  in  the  energy  bin  I'  produces  a  count 
C  in  the  pulse-height  bin  J'. 

C 

DO  170  L  =  1,NRESP,3 
J  =  J  +  1 

DO  160  K  =  1,NBINS 

FINRESP(J,K)  =  (SUBRESP(L,K)  SUBRESP(L+1,K) 

1  +  SUBRESP(L+2,K))/3.0 

160  CONTINUE 
170  CONTINUE 
C 

C  Output  the  response  matrix  to  a  file  for  ease  of  reading— BGRESP. OUT. 
C 

OPEN(  UNIT=2 ,  FILE= '  BGRESP .  OUT ' ,  STATUS^:  ■  , ' ) 

WRITE(2,180) 

180  FORHAT(T10,'**ADOPTED  BINS  FOR  THE  RESPONSE  HATRIX**'/) 
WRITE(2,190) 

190  FORHAT(lX,54('-')/lX,54('-')) 

WRITE(2,200) 

200  FORHAT(T6,'Bin',T15/Bin',T27,'Bin',T44/Bin', 

1  /T5, 'Number' .T13, 'Centroid', T26, 'Width', T43, 

2  'Limits', /il^,'(keV)',T26,'(keV)',T43,'(keV)') 

WRITE(2,190) 

DO  220  I  =  1,NBINS 

BINC  =  (PHLIHS(I)  +  PHLIHS(I+l))/2.0 
BINW  =  PHLIHS(I+1)  -  PHLIHS(I) 

WRITE(2,210)  I,BINC,BINW,PHLIHS(I),PHLIHS(I+1) 

210  FORHAT(T6,I3,T13,F8.2,T25,F7.2,T36,F8.2, 

1  1X,'-',1X,F8.2) 

220  CONTINUE 
WRITE(2,190) 

WRITE (2,230)  TYPE 

230  FORHAT(//T10, 'Response  Hatrix  For  ',A17,'  BGO  Inorganic', 

1  '  Scintillator'/T19, 'Column  =  Gamma-ray  Energy,  Row  = 

2  'Pulse  Height') 

WRITE(2,240) 
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240  FOmT(//) 

DO  270  I  =  1,NBINS,7 
H«  =  NBINS  -  I 
IF(HM.GT.6)  THEN 
NN  =  6 

ELSE 

NN  =  HH 

ENDIF 

DO  260  J  =  1, NBINS 

WliITE(2,250)  (FINKESP(I+K,J),K=0,NH) 

250  F0RMAT(1X,7(2X,1PE9.3)) 

260  CONTINUE 

WEITE(2,240) 

270  CONTINUE 

CLOSE(UNIT=2) 

C 

C  Output  the  file  in  a  succint  fornat  for  input  to  an  unfolding  code. 

C  Note  that  the  highest  energy  response  of  the  natrix,  I=NBINS,  is  not 
C  printed  to  this  file.  The  unfolding  code  intended  for  use  of  this 
C  file,  BGO.HAT,  does  not  need  the  highest  energy  response  of  the 
C  natrix.  If  a  different  unfolding  code  is  used,  nodification  of 
C  this  file  nust  be  nade. 

C 

0PEN( UNIT=2 , FILE= 'BGO. NAT ' , STATUS^ 'NEft ' ) 

WRITE(2,280)  ((FINPESP(I,J), 1=1, NBINS-1),J=1, NBINS) 

280  F0RHAT(8(1X,1PE9.3)) 

CL0SE(UNIT=2) 

C 

C  Output  the  number  of  bins  in  the  response  matrix  and  the  pulse- 
C  height  values  of  these  bins  to  the  file  BGLIHS.DAT.  This  file 
C  is  used  by  an  unfolding  code.  The  pulse-height  values  in  this 
C  file  are  in  HeV  and  not  keV. 

C 

0PEN(UNIT=2,FILE='BGLIHS.DAT',STAT0S='NEW') 

WRITE(2,290)  NBINS 
290  F0RHAT(I5) 

WkITE( 2, 300)  ( PHLIHS( I ) /lOOO .0,1=1 ,NBINS+1 ) 

300  FORNAT(5(6X,F7.4)) 

CLOSE(UNIT=2) 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  BGIRESP(ENRESP) 

C 

C  Subroutine  BGIRESP  calculates  an  absolute  differential  efficiency 
C  curve  for  gamma  rays  for  a  2.54-cn  x  2.54-cn  BGO  inorganic  scintillator 
C  for  any  specified  gamma  ray  energy.  This  curve  is  what  is  normally 
C  called  the  response  function.  The  response  function  is  generated 
C  from  an  interpolation  scheme  using  12  data  files  which  contain 
C  experimental  response  functions  for  monoenergetic  fluxes  of  gamma 
C  rays  incident  on  the  detector  (source-to-detector  distance  eguals 
C  one  meter)  with  the  following  energies:  123.6,  279.2,  514.0,  661.6, 

C  834.8,  1115.5,  1274.5,  1836.0,  2754.0,  4440.0,  6140.0,  and 
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C  11670.0  keV.  The  calculated  response  function  is  stored  in  the  array 
C  RESP.  RESP(r)  is  the  probability  per  unit  pulse  height  (in  keV) 

C  for  one  ganna  ray  of  the  specified  energy  to  produce  a  count  in  the 
C  pulse-height  bin  I'. 

C 

C  Subroutine  Suimary: 

C 

C  Other  Subroutines  Used: 

C 

C  None. 

C 

C  Functions  Used: 

C 

C  None. 

C 

C  CoBBOn  Blocks  Used: 

C 

C  ARRAYS 

C 
C 
C 

C  DATARAYS 

C 
C 

C  PARK 

C 
C 
C 

C  Local  Variables  Used: 


==  Contains  8  1-D  arrays  shared  between  this 
subroutine  and  the  subroutine  BG3RESP  which 
is  identical  to  this  subroutine  except  for  the 
data  statenents. 

==  Contains  2  1-D  arrays  shared  between  this 
subroutine  and  the  Bain  prograB  and  subroutine 
BG3RESP. 

==  Contains  3  parameters  shared  between  this 
subroutine  and  the  main  program  and  subroutine 
BG3RESP. 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


BGNORH  ==  Real  1-D  12  element  array.  Contains  the 

normalization  i  for  each  experimental  response 
function  data  file  in  order  to  convert  the  data 
from  relative  counts  to  counts  per  incident 
gamma  ray. 

CTSHI  =-  Real  l-D  1024  element  array.  In  Common  Block 
ARRAYS.  Same  as  array  CTSLO  except  that  it 
contains  the  information  for  the  gaima-ray 
energy  which  is  the  closest,  but  greater  than 
or  equal  to,  uie  'sVecified  energy. 

CTSLO  ==  Real  1-D  1024  element  array.  In  Common  Block 
ARRAYS.  Contains  the  experimental  response 
function  (pulse-height  integrated— counts/ 
gaBBa  ray)  for  the  qamna-ray  energy  which 
is  the  closest,  but  less  than  or  equal  to, 
the  energy  passed  to  the  subroutine  for 
response  function  computation. 

DELTAfll  ==  Real  variable.  The  absolute  energy  distance 
between  EHI  and  ENRESP. 

DELTALO  ==  Real  variable.  The  absolute  energy  distance 
between  ELO  and  ENRESP. 

WTOEHI  ==  Real  1-D  1024  element  array.  In  CoBBon  Block 
ARRAYS.  CTSHI  divided  by  the  experinental 
bin  width  and  multiplied  by  BGNORN  so  that  the 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


experinental  data  is  counts/keV/gamiia  ray. 

DNDELO  ==  Real  1-D  1024  elenent  array.  In  Comaon  Block 
ARRAYS.  CTSLO  divided  by  the  experinental 
bin  width  and  nultiplied  by  BGNORM. 

DUHVAL  ==  Real  variable.  A  dunmy  variable. 

EGAKKA  ==  Real  1-D  12  elenent  array.  Contains  the  ganna-ray 

energy  values  for  the  12  neasured  response  functions. 

EHI  ==  Real  variable.  The  energy  of  the  data  file  which 
brackets  ENRESP  on  the  high  side. 

ELO  ==  Real  variable.  The  energy  of  the  data  file  which 
brackets  ENRESP  on  the  low  side. 

ENRESP  ==  Real  variable.  Passed  to  the  subroutine  through 
the  call.  This  is  the  gaana-ray  energy  for 
conputation  of  a  response  function. 

GFILE  ==  Character  1-D  12  elenent  array.  Contains  the 

file  nanes  of  the  experinental  response  function 
data  files. 

GSLOPE  ==  Real  1-D  12  elenent  array.  Contains  the  slope 
in  keV/Channel  f  for  the  pulse-height  axis 
calibration  for  each  experinental  response 
function  data  file. 

GYINT  ==  Real  1-D  12  elenent  array.  Contains  the  intercept 
in  keV  for  the  pulse-height  axis  calibration 
for  each  experinental  response  function  data  file. 

INC  ==  Integer  variable.  A  counter. 

JHIl  ==  Integer  variable.  A  counter. 

JHI2  ==  Integer  variable.  A  counter. 

JLOl  ==  Integer  variable.  A  counter. 

JL02  ==  Integer  variable.  A  counter. 

HBI  ==  Integer  variable.  A  counter. 

HLO  ==  Integer  variable.  A  counter. 

PHBIN  ==  Real  variable.  In  Connon  Block  PARH.  See 
definition  in  the  nain  progran. 

PHCUT  ==  Real  1-D  12  element  array.  One  value  for  each 
experinental  response  function  data  file.  Above 
this  pulse-height  value,  the  data  file  is  shifted 
downward  or  upward  in  pulse  height  in  order  to 
match  the  pulse  height  of  the  response  function 
for  the  gamma-ray  energy  ENRESP. 

PHVAL  ==  Real  variable.  A  dummy  variable  for  storing 
pulse-height  values. 

PINCHI  ==  Real  variable.  The  pulse-height  increments  needed 
to  span  the  pulse-height  range  WIDTHHI. 

PINCLO  ==  Real  variable.  The  pulse-height  increments  needed 
to  span  the  pulse-height  range  WIDTHLO. 

PLCUT  ==  Real  1-D  12  element  array.  One  value  for  each 
experimental  response  function  data  file.  Below 
this  pulse-height  value,  the  pulse-height  axis 
of  the  data  file  is  not  modified. 

PRESPH  ==  Real  variable.  In  Common  Block  PARN.  Same  as 
PRANE  in  the  main  program. 

PRESPL  ==  Real  variable.  In  Common  Block  PARH.  Same  as 
PRANL  in  the  main  program. 

PVALHI  ==  Real  1-D  1024  element  array.  In  Common  Block 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 


c 

c 

c 

c 
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ARRAYS.  Contains  the  pulse-heiqht  bin  values 
(in  keV)  used  to  generate  the  response  function 
stored  in  RESPHI. 

PVALLO  ==  Real  1-D  1024  element  array.  In  CoBBon  Block 
ARRAYS.  Contains  the  pulse-heiqht  bin  values 
(in  keV)  used  to  generate  the  response  function 
stored  in  RESPLO. 

RESP  ==  Real  1-D  1000  eleaent  array.  In  CoMion  Block 
DATARAYS.  See  definition  in  the  main  program. 

RESPHI  ==  Real  1-D  1000  element  array.  In  Common  Block 
ARRAYS.  Contains  the  same  response  function 
as  stored  in  DNDEHI  except  it  has  been  shifted 
and  compressed  along  the  pulse-heiqht  axis 
in  order  to  simulate  the  pulse-hei^t  axis 
of  the  response  function  for  the  gamma-ray  energy 
ENRESP. 

RESPLO  ==  Real  1-D  1000  element  array.  In  Common  Block 
ARRAYS.  Contains  the  same  response  function 
as  stored  in  DNDELO  except  it  has  been  shifted 
and  stretched  along  the  pulse-heiqht  axis  in 
order  to  simulate  the  pulse-heiqht  axis  of  the 
response  function  for  the  gamma-ray  energy  ENRESP. 

RESPPH  ==  Real  1-D  1000  element  array.  In  Common  Block 
DAT.ARAYS.  See  definition  in  the  main  program. 

WIDTHFI  ==  Real  variable.  Range  of  data  in  the  experimental 

respc.'.se  function  data  file  for  the  gamma-ray  energy 
EHI  which  must  be  compressed. 

WIDTHLO  ==  Real  variable.  Range  of  data  in  the  experimental 

response  function  data  file  for  the  gamma-ray  energy 
ELO  which  must  be  stretched. 

Set  up  the  common  blocks  used  between  this  subroutine  and  the  main 
program  and  the  subroutine  BG3RESP. 

COfflON/PARM/PRESPL , PRESPH , PHBIN 

COHHOH/DATARAYS/RESP( 1000) ,RESPPH( 1000) 

COHNON/ARRAYS/CTSLO( 1024 ) ,CTSHI ( 1024 ) ,DNDEL0(1024 ) , 

1  DNDEHI ( 1024 ) , PVALLO) 1024 ) , PVALHI ( 1024 ) , 

2  RESPLO(  1000  i,R13Pn(  1000) 

Dimension  the  variables. 

CHARACTER  GFILE(12)*11 

DIMENSION  GSLOPE(12),GYINT(12),BGNORM(12),PLCUT(12), 

1  PHCUT(12),EGAMMA(12) 

Load  the  12-element  arrays  with  the  necessary  data.  GSLOPE,  GYINT, 
PLCUT,  PHCOT,  and  EGAHNA  are  energy  or  pulse-height  values  and 
the  numbers  all  have  units  of  keV. 

DATA  GFILE/ ' BG124 . DIT ' , ' BG279 . DIT ' , ' BG5i4 . DIT' , 'BG662 . DIT' , 

1  'BG835.D1T' , 'BG1116.D1T' , 'BG1275.D1T' , 

2  'BG1836.D1T' , 'BG2754.D1T' , 'BG4440.D1T' , 

3  'BG6140.D1T','BG11670.D1T7 
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DATA  GSLOPE/0. 41803, 0.42656, 0.69732, 0.8445, 1.05882, 1.4503, 

1  1.68546,2.2352,3. 17087 , 4 . 90506,7. 99978 , 13 . 75606 / 

DATA  GYINT/11. 1417, 12. 0048, 10. 55345,8. 559,8. 7637,6. 5508, 

1  12.45844,10.786,-6.14645,-3.7977,-23.565,-23.132/ 

DATA  BGNORM/2 . 26E-06 , 1 . 1 4E-05, 3 . 096E-06 , 7 . 16E-07 , 1 . 49E-06 , 

1  1 . 87E-06 , 1 . OlE-06 , 2 . 13E-06 , 1 . 59E-06, 5 . 62E-07 , 

2  8.94E-07,1.48E-05/ 

DATA  PLCUT/67. 0,107. 0,243. 0,255.0, 300.0,317.0, 330. 0,660.0, 

1  700.0,733.0,650.0,700.0/ 

DATA  PHCUT/90. 0,200. 0,342. 0,425. 0,573. 0,803.0,960.0, 

1  1158.0,1560.0,3167.0,4860.0,9583.0/ 

DATA  EGAHMA/123.6,279.2,514.0,661.6,834.8,1115.5,1274.5, 

1  1836.0,2754.0,4440.0,6140.0,11670.0/ 

C 

C  This  loop  determines  which  experimental  response  function  data  files 
C  bracket  (in  enerqy)  the  response  function  which  is  to  be  computed 
C  at  the  gamma-ray  energy  ENBESP. 

C 

DO  10  I  =  2,12 

IF(ENPESP.LE.EGAMHA(I))  THEN 
ELO  =  EGAfflA{I-l) 

EHI  =  EGAHHA(I) 

HLO  =  I  -  1 
Hfil  =  I 
GO  TO  15 

ENDIP 

10  CONTINUE 
15  CONTINUE 
C 

C  Open  the  two  experimental  response  function  data  files  which  bracket 
C  the  response  function  to  be  calculated  and  read  the  data— the  data 
C  is  such  that  it  has  the  units  of  counts.  That  is,  each  experimental 
C  response  function  data  file  is  pulse-height  integrated  (counts  in  a 
C  bin)  and  is  not  normalized  to  the  total  number  of  gamma-rays  incident 
C  on  the  detector  to  produce  the  data  file. 

C 

OPEN(  UNIT=2 ,  FILE=GFILE  ( HLO ) ,  STATUS= 'OLD' ) 
0PEN(0NIT=3,FILE=GFILE(HHI) ,STATUS='OLD' ) 

EEAD(2,20)  (CTSLO(I),I  =  1,1024) 

20  F0RMAT(1X,6F13.2) 

READ(3,20)  (CTSHI(I),I  =  1,1024) 

CLOSE (UNIT=2) 

CLOSE(UNIT=3) 

C 

C  Normalize  the  experimental  response  function  and  diviu  jh  binned 
C  count  by  the  pulse-height  bin  width.  Thus,  DNDELO  and  DNDEHI  are 
C  experimental  absolute  differential  efficiency  curves— i.e.,  they 
C  have  units  of  counts/unit  pulse  height  (keV) /incident  gamma  ray. 

C 

DO  30  I  =  1,1024 

DNDELO(I)  =  CTSLO(I)*BGNORH(KLO)/GSLOPE(HLO) 

DNDEHI(I)  =  CTSHI(I)*BGNORH(HHI)/GSLOPE(HHI) 

30  CONTINUE 
C 


59 


C  Determine  the  difference  in  energies  between  the  response  function 
C  to  be  calculated  and  the  two  response  functions  bracketing  this 
C  response  function  and  which  will  be  used  for  the  interpolation 
C  process. 

C 

DELTALO  =  ENRESP  -  ELO 
DELTAHI  =  EHI  -  ENRESP 
C 

C  The  interpolation  process  is  as  follows:  first,  the  experimental 
C  response  function  which  is  for  the  gamma-ray  energy  just  below  or 
C  or  equal  to  ENRESP  is  modified.  This  modification  is  a  modification 
C  of  the  pulse-height  axis  only.  The  pulse-height  axis  is  the  original 
C  pulse-height  axis  below  the  value  PLCUT,  shifted  upward  by  the  amount 
C  DELTALO  a^ve  the  value  PHCUT,  and  stretched  in  between  to  fill 
C  the  remainder  of  the  pulse-height  axis.  A  similar  modification 
C  is  made  to  the  experimental  response  function  which  is  for  the 
C  gamma-ray  energy  just  above  or  equal  to  ENRESP—the  pulse-height  axis 
C  is  the  original  below  the  value  PLCUT,  shifted  downward  by  the  amount 
C  DELTAHI  above  the  value  PHCUT,  and  compressed  in  between  to  fill 
C  the  remainder  of  the  pulse-height  axis.  Second,  these  two  pulse- 
C  height  modified  response  functions  are  then  rebinned  into  1000 
C  data  points  wnrcii  are  equally  spaced  about  the  pulse-height  range 
C  specified  by  the  user  for  the  response  matrix,  i.e.,  the  range 
C  PRESPL  to  PRESPH.  Third,  these  resulting  two  experimental  response 
C  functions  are  then  used  to  generate  a  response  function  at  the 
C  gamma-ray  energy  ENRESP  through  a  log-log  interpolation  along 
C  the  y  axis  ( counts /keV/qaiia  ray  axis). 

C 

INC  =  0 
C 

C  First  step  of  the  interpolation  process. 

C 

DO  40  I  =  1,1024 

DUHVAL  =  GSLOPE(HLO)*I  +  GYINT(HLO) 

IF(DOHVAL.LE.PLCUT(MLO))  THEN 
PVALLO(I)  =  DUHVAL 
JLOl  =  I 

ENDIF 

IFiDUfVAL.GE.PHCUKHLO)  T  'H 
INC  =  INC  +  1 

PVALLO(I)  =  DUHV’AL  +  DELTALO 

ENDIF 

40  CONTINUE 

JL02  =  1024  -  INC  +  1 

WIDTHLO  =  PVALL0(JL02)  -  PVALLO(JLOi) 

PINCLO  =  WIDTHLO/ (JL02  -  JLOl) 

DO  50  I  =  ju01+l,JL02-l 

PVALLO(I)  =  PVALLO(I-l)  +  PINCLO 
50  CONTINUE 
INC  =  0 

DO  60  I  =  1,1024 

DUHVAL  -  GSLOPE(HHI)*I  +  GYINT(HHI) 

IF(DUHVAL.LE.PLCUT(HHI))  THEN 
PVALHKI)  =  DUHVAL 
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JHIl  =  I 

ENDIF 

IF(DUHVAL.GE.PHCUT(HHI))  THEN 
INC  =  INC  +  1 

PVALHI(I)  =  DUHVAL  -  DELTAHI 

ENDIF 

60  CONTINUE 

JHI2  =  1024  -  INC  +  1 

WIDTHHI  =  PVALHI(JHI2)  -  PVALHI(JHIl) 

PINCHI  =  WIDTHHI /(JHI2  -  JHIl) 

DO  70  I  =  JHI1+1,JHI2-1 

PVALHI(I)  =  PVALHI(I-l)  +  PINCHI 
70  CONTINUE 
C 

C  Second  step  of  the  interpolation  process. 

C 

PHBIN  =  (P8ESPH  -  PRESPL)/1000.0 
PHVAL  =  PRESPL  -  PHBIN/2.0 
DO  100  I  =  1,1000 

PHVAL  =  PHVAL  +  PHBIN 
DO  80  J  =  2,1024 

IF(PHVAL.LE.PVALLO(J))  THEN 

RESPLO(I)  :=  (DNDELO(J)  -  DNDELO(J-l))* 

1  (PHVAL  -  PVALLO(J-l))/ 

2  (PVALLO(J)  -  PVALLO(J-l))  + 

3  DNDELO(J-l) 

IF(RESPLO(I).LT.C.O)  EESPLO(I)  =  0.0 
GO  TO  85 

ENDIF 

80  CONTINUE 

85  CONTINUE 

DO  90  J  =  2,1024 

IF(PHVAL.LE.PVALHI(J))  THEN 

RESPHI(I)  =  (DNDEHI(J)  -  DNDEHI(J-l))* 

1  (PHVAL  -  PVALHI(J-l))/ 

2  (PVALHI(J)  -  PVALHI(J-l))  t 

3  DNDEHKJ-l) 

IF(RESPHI(I).LT.0.0)  RESPHI(I)  =  0.0 
GO  TO  95 

ENDIF 

90  CONTINUE 

95  CONTINUE 

100  CONTINUE 
C 

C  Third  and  final  step  of  the  interpolation  process. 

C 

DO  110  I  =  1,1000 

IF(RESPLO(I).EQ.O.O.OR.RESPHI(I).EQ.O.O)  THEN 
RESP(I)  =  RESPLO(I)*(LOG(EHI)  -  LOG(ENRESP)) 

1  /(LOG(EHI)  -  LOG(ELO))  + 

2  RESPHI(I)*(LOG(ENRESP)  -  LOG(ELO)) 

3  /(LOG(EHI)  -  LOG(ELO)) 

ELSE 

RESP(I)  =  LOG(RESPLO(I))*(LOG(EHI)  -  LOG(ENRESP)) 
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1  ,'(LOG(EHI)  -  LOG(ELO))  + 

2  LOG(RESPHI|I))*(LOG(EHRESP)  -  LOG(ELO)) 

3  /(LOG(EHI)  -  LOG(ELO)) 

RESP(I)  =  EXP(RESP(I)) 

ENDIF 

RESPPH(I)  =  PHBIN*I  t  PRESPL  -  PHBIN/R.O 
no  CONTINUE 
RETURN 
END 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCcCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  BG3RESF(ENRESP) 

C 

C  Subroutine  BG3RESe  calculates  an  absolute  ditterential  efficiency 
C  curve  for  qaima  rays  for  a  7.62-cn  x  7.62-co  BGO  inorganic  scintillator 
C  for  any  specified  gama  ray  energy.  This  subroutine  is  identical 
C  to  BGIRESP.  See  the  subroutine  BGIRESP  for  docunentation. 

C 

COHHON/''\EN/PRESPL,PRESPH,PHBIN 
COmON/DATARAYS/RESPi  1000)  ,RESPPH(  1000) 

COHNON/ARRAYS/CTSIO( 1024 ) ,CTSHI { 1024 ) ,DNDELO( 1024 ) , 

1  DNDEfll ( 1024 ) , PVALLO) 1024 ) , PVALHI ( 1024 ) , 

2  RESPLO) 1000 ),RESPHI( 1000) 

CHARACTER  GFILE(12)*11 

DIMENSION  GSLOPE(12),GYINT(12),BGNORM(12),PLCUT(12), 

1  PHCUT(12),EGAHMA(12) 

DATA  GFILE/'BG124.D3T'/BG279.D3T','BG514.D3T','BG662.D3T', 

1  'BG835.D3T' , 'BG1116.D3T' , 'BG1275.D3T' , 

2  'BG1836.D3T' , 'BG2754 .D3I' , 'BG4440.D3T' , 

3  'BG6140.D3I'/BG11670.D3T'/ 

DATA  GSLOPE/0. 40024, 0.3983, 0.66238, 0.90496, 1.03376, 1.4835, 

1  1.67676,2.34348,3.1788,5.31817,7.50905,13.29843/ 

DATA  GYINT/12. 7296, 9. 7776, 7. 9805, 5. 43216, 4. 5459,0.11605, 

1  -6.0159,-17.099,-19.1988,-13.3556,-31.3755, 

2  -17.4087/ 

DATA  BGNORH/2 . 51E-07 , 1 . 92E-06 , 3 . 4 4 E-07 , 7 . 96E-08 , 1 . 65E-07 , 

1  2.08E-07,1.12E-07,2.36E-07,5.44E-07,1.66E-07, 

2  1.37E-07,2.55E-0t 

DATA  PLCUT/  7  0, 10’. 0,260.0, 25/. C, 330.0, 3 - J.0, 340.0,601.0, 

1  601.0,750.0,700.0,700.0/ 

DATA  PHCUT/90. 0,200. 0,357. 0,470. 0,605. 0,880. 0,612.0, 

1  1129.0,1562.0,3233.0,4900.0,9370.0/ 

DATA  EGAMHA/123.6,279.2,514.0,661.6,834.8,1115.5,1274.5, 

1  1836.0,2754.0,4440.0,6140.0,11670.0/ 

DO  1C  I  =  2,12 

IF{ENRESP.LE.EGAHHA(I))  THEN 
ELO  =  EGAHHA(I-l) 

EHI  =  EGAHHA(I) 

HLO  =  I  -  1 
HHI  =  I 
GO  TO  15 

ENDIF 

10  CONTINUE 
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15  CONTINUE 

OPEN ( UHIT-2 , riEE=GriLE ( HLO ) , STATUS= 'OLD' ) 
OPEN(UNIT=3,FILE=GFILE(HHI),STATUS='OLD') 

READ! 2, 20)  (CTSLO(I),I  =  1,1024) 

20  FORMAT(lX,6Fi3.2) 

EEAD(3,20)  (CT?Hr(I),I  =  1,1024) 

CLOSE (UNIT=2) 

CLOSE(UNIT=3) 

DO  30  I  ^  1,1024 

DNDELO(I)  =  CTSLO(I)*BGNORH(HLO)/GSLOPE(MLO) 
DNDEHI(I)  =  aSHI(I)*BGNORM(HHI)'GSLOPE(NHI) 
30  CONTINUE 

DELTALO  =  ENRESP  -  ELO 
DELTAHI  =  EHI  -  ENRESP 
INC  =  0 

DO  40  I  =  1,1024 

DUIf/AL  =  GSLOPE(HLO)*I  +  GYINT(HLO) 
IF(DUMVAL.LE.PLCUT(HLO))  THEN 
PVALLO(I)  =  DUHV’AL 
JLOl  =  I 

ENDIF 

IF(DU«VAL.GE.PHCUT(HLO))  THEN 
INC  =  INC  +  1 

PVALLO(I)  =  DUHVAL  +  DELTALO 

ENDIF 

40  CONTINUE 

JL02  =  1024  -  INC  +  1 

HIDTHLO  =  PVALL0(JL02)  -  PVALLO(JLOl) 

PINCLO  =  WIDTHL0/(JL02  -  JLOl) 

DO  50  I  =  JL01+1,JL02-1 

PVALLO(I)  =  PVALLO(M)  +  PINCLO 
50  CONTINUE 

INC  =  0 

DO  60  I  =  1,1024 

DUHVAL  ^  GSLOPE(HHI)*I  +  GYINT(HHI) 

IF( DUHVAL. LE.PLCUT(HHI))  THEN 
PVALHI(I)  =  DUHVAL 
JHIl  ^  I 

ENDIF 

IF(DUHVAL.GE.PflCUT(HHD)  THEN 
INC  =  INC  M 

PVALHI(I)  =  DUfT.’AL  -  DELTAHI 

ENDIF 

60  CONTINUE 

JHI2  =  1024  -  INC  T  1 

WIDTHHI  =  PVALHI(JHI2)  -  PVALHI(JHIl) 

PINCHI  =  WIDTHHI/(JHI2  -  JHIl) 

DO  70  I  =  JHI1+1,JHI2-1 

PVALHI(I)  =  PVALHI(I-l)  *  PINCHI 
70  CONTINUE 

PHBIN  =  (PRESPH  -  PRESPL)/ 1000.0 
PHVAL  ^  PRESPL  -  PHBIN/2.0 
DO  100  I  =  1,1000 

PffVAL  ^  PPAL  +  PHBIN 
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DO  80  J  =  2,1024 

if(phval.le.pvallo(J))  then 

RESPLO(I)  =  (DNDELO(J)  -  DNDELO(J-l) )* 

1  (PHVAL  -  PVALLO(J-l))/ 

2  (PVALLO(J)  -  PVALLO(J-l))  + 

3  DNDELOlJ-1) 
IF(RESPLO(I).LT.O.O)  RESPLOd)  =  0.0 
GO  TO  85 


80 

85 


90 

95 

100 


no 


ENDIF 
CONTINUE 
CONTINUE 
DO  90  J  =  2,1024 

IF(PHVAL.LE.PVALHI(J)|  THEN 

RESPHI(I)  -  (DNDEHKJ)  -  DNDEUI(J-l) )* 

1  (PHVAL  -  PVALHI(J-l))/ 

2  (PVALHKJ)  -  PVALHI(J-l))  + 

3  DNDEHKJ-l) 

IF(RESPHI(I).1J.0.0)  RESPHI(I)  =0.0 
GO  TO  95 

ENDIF 

CONTINUE 

CONTINUE 

CONTINUE 

DO  no  I  =  1,1000 

IF(RESPLO(I).EQ.O.O.OR.RESPHI(I).EQ.O.O)  THEN 
RESP(I)  =  RESPLO|I)*(LOG(EHI)  -  LOG(ENRESP)) 

1  '(LOG(EHIi  -  LOG(ELO))  + 

2  RESPHI(I)*(LOG(ENRESP)  -  LO€(ELO)) 

3  ,  (  LOGiEHD  -  LOGIELO)) 

ELSE 

RESP(I)  =  LOG(RESPLO(I)l*(LOG(EHI)  -  LOG(ENRESP)) 

1  ;(LOG(EHI)  -  LOG(ELO))  + 

2  LOG(RESPHI(I)l*(LOG{ENRESP)  -  LOG(ELO)) 

3  /(LOG(EHI)  -  LOG(ELO)) 

RESP(I)  =  EXP(RESP(I)) 

ENDIF 

RESPPH(I)  =  PHBIN*I  ^  PRESPL  -  PHBIN/2.0 
CONTINUE 
RETURN 
END 


C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

FUNCTION  TERP ( PHVAL , PHITE , COUNTS  t 
C 

C  This  function  is  a  simple  interpolation  routine  and  is  self-explanatory. 

C 

DIMENSION  PHITEI 1000 1 , COUNTS! 1000 ) 

IfiPHVAL.LE.PHITEdl)  THEN 
TERP  =  COUNTS! 1 1 
RETURN 


ENDIF 

IF(PHVAL.GT.PHITE(1000|)  THEN 
TERP  =  COUNTS! 1000) 
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EETURN 


ENDIF 

DO  10  I  =  2,1000 

IF(PHVAL.LE.PHITE(I))  THEN 

TERP  =  (PPAL  -  PHITE(I-1))/(PHITE(I)  -  PHITE(I-l)) 
1  *(COUNTS(I)  -  COUHTS(I-l))  +  COUNTS(I-l) 

RETURN 

ENDIF 

10  CONTINUE 
END 
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